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ABSTRACT 

The combination of gravitational lensing and stellar dynamics breaks the mass-anisotropy degeneracy 
and provides stringent constraints on the distribution of luminous and dark matter in early- type (E/SO) 
galaxies out to z w 1. We present new observations and models of three lens systems (CFRS03.1077, 
HST14176+5226, HST15433+5352) and the combined results from the five field E/SO lens galaxies at 
z as 0.5—1.0 analyzed as part of the Lenses Structure & Dynamics (LSD) Survey. Our main results 
are: (i) Constant mass-to-light ratio models are ruled out at > 99% CL for all five E/SO galaxies, 
consistent with the presence of massive and extended dark-matter halos. The range of projected dark- 
matter mass fractions inside the Einstein radius is /dm=0.37-0.72, or 0.15-0.65 inside the effective 
radius R c for isotropic models, (ii) The average effective power-law slope of the total (luminous+dark; 
Ptot oc r~ 7 ) mass distribution is (7') = 1.75±0.10 (1.57±0.16) for Osipkov-Merritt models with anisotropy 
radius ri — oo(R e ) with an rms scatter of 0.2 (0.35), i.e. marginally flatter than isothermal (7' = 2). 
The ratio between the observed central stellar velocity dispersion and that from the best-fit singular 
isothermal ellipsoid (SIE) lens model is (/sie) = (c/csie) = 0.87 ± 0.04 with 0.08 rms, consistent with 
flatter-than-isothermal density profiles. Considering that 7' > 2 and /sie > 1 have been reported for 
other systems (i.e. B1608+656 and PG1115+080), we conclude that there is a significant intrinsic scatter 
in the slope of the mass-density profile of lens galaxies (rms ~ 15%) , similar to what is found for local E/SO 
galaxies. Hence, the isothermal approximation is not sufficiently accurate for applications that depend 
critically on the slope of the mass-density profile, such as the measurement of the Hubble Constant 
from time-delays, (iii) The average inner power-law slope 7 of the dark-matter halo is constrained to 
be (7) = 1-31q 4 (68% C.L.), if the stellar velocity ellipsoid is isotropic (r, = 00) or an upper limit of 
7 < 0.6, if the galaxies are radially anisotropic (r, = i?, c ). The observed range of slopes of the inner 
dark-matter distribution is consistent with the results from numerical simulations only for an isotropic 
velocity ellipsoid and if baryonic collapse and star- formation do not steepen dark-matter density profiles, 
(iv) The average stellar mass-to-light ratio evolves as d\og{M*/L-a)/dz=— 0.72 ± 0.10, obtained via a 
Fundamental Plane analysis. An independent analysis based on lensing and dynamics gives an average 
(d\og(M^/L-B)/dz)—— 0.75 ± 0.17. Both values indicate that the mass-to-light ratio evolution for our 
sample of field E/SO galaxies is slighly faster than those in clusters, consistent with the hypothesis that 
field E/SO galaxies experience secondary bursts (~10% in mass) of star formation at z < 1. These 
findings are consistent with pure luminosity evolution of E/SO galaxies in the past 8 Gyrs, and would be 
hard to reconcile with scenarios involving significant structural and dynamical evolution. 

Subject headings: gravitational lensing — galaxies: elliptical and lenticular, cD — galaxies: evolution 
— galaxies: formation — galaxies: structure 



1. INTRODUCTION 

A central assumption of the current standard ACDM 
cosmological model is that galaxies form and evolve inside 
dark-matter halos (White & Rees 1978; Blumenthal et al. 
1984; Davis et al. 1985). Dark matter halos are ubiquitous 
and possibly universal in their density profile (Navarro, 
Frenk & White 1996,1997, hereafter NFW; Moore et al. 
1998) and they dominate the dynamics of large scale struc- 
tures. In spite of decades of searches and technological 
advances our empirical knowledge of dark halos remains 



very sparse. 

In the local Universe, dark-matter halos have been con- 
vincingly detected - predominantly through dynamical 
tracers - in spiral galaxies (e.g. Rubin, Thonnard, & Ford 
1978,1980; Faber & Gallagher 1979; van Albada & Sancisi 
1986; Salucci & Burkert 2000; Jimenez, Verde & Oh 2003), 
dwarf and low surface brightness galaxies (de Blok & Mc- 
Gaugh 1997; van den Bosch et al. 2000; Swaters, Madore 
& Trewella 2000), clusters of galaxies (Zwicky 1937; Kneib 
et al. 1993; Lombardi et al. 2000; Sand et al. 2002, 2004; 
Kneib et al. 2003), and - at least in some cases - in early- 
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type type galaxies (e.g. Fabbiano 1989; Mould et al. 1990; 
Matsushita et al. 1998; Loewenstein & White 1999; Saglia, 
Bertin & Stiavelli 1992; Bertin et al. 1994; Arnaboldi et 
al. 1996; Franx et al. 1994; Carollo et al. 1995; Rix et 
al. 1997; Gerhard et al. 2001; Borricllo, Salucci & Danese 
2003; Seljak 2002). However, observational evidence re- 
garding dark-matter halos of early- type galaxies (E/SOs) 
is limited. In a number of cases, constant mass-to-light 
models, with a mass-to-light ratio consistent with those of 
a normal stellar population, appear sufficient to explain 
the information available from mass tracers, and there is 
no need to invoke the existence of dark-matter halos (e.g. 
Bertin et al. 1994; Romanowsky et al. 2003). 

If detection of dark-matter is hard and often ambigu- 
ous, decomposing the mass distribution into a luminous 
(mostly stellar) and dark-matter component, to measure 
their relative contribution and spatially resolved proper- 
ties, has been possible in only very few cases with varying 
results (see references above). The main hurdles to over- 
come are the paucity of dynamical tracers at large radii 
(such as HI gas in spirals) and the degeneracy between 
kinematic properties of dynamical tracers (e.g. anisotropy 
for stellar dynamics) and the mass distribution. For sim- 
plicity, we will refer to the latter problem as the mass- 
anisotropy degeneracy. 

The distant Universe (z > 0.1) is an almost completely 
uncharted territory. Gravitational lensing has provided 
evidence for a mass distribution that is more extended 
than the stellar component, either by the analysis of in- 
dividual systems or by considering statistical ensembles 
(e.g. Kochanek 1995, Rusin & Ma 2001, Ma 2003, Rusin, 
Kochanek & Keeton 2003; Cohn et al. 2001, Munoz, 
Kochanek & Keeton 2001; Rusin et al. 2002; Winn, Rusin 
& Kochanek 2003; Wucknitz, Biggs & Browne 2004). Un- 
fortunately, for most lenses it has proven very difficult to 
reliably separate the luminous from the dark-matter com- 
ponent using lensing alone, and determine a precise dark- 
matter density profile and mass fraction. 

In spite of the difficulties, a detailed exploration of high- 
z galaxies would come with a great reward, offering the op- 
portunity to map directly the evolution of dark and stellar 
mass over cosmic time. By mapping the time (i.e. redshift) 
evolution of the relative distributions of luminous and dark 
matter in early-type galaxies - as well as the evolution of 
stellar mass-to-light ratio and the slope of dark-matter ha- 
los - we can address directly the following questions. How 
and when is mass assembled to form early-type galaxies? 
What is their star formation history? Are dark-matter ha- 
los characterized by cuspy mass profiles in the center as 
predicted by numerical simulations? What is the role of 
star formation in shaping the total mass distribution of 
early-type galaxies? Do isolated early-type galaxies un- 
dergo internal structural/dynamical evolution? 

To answer these questions, we are undertaking the 
Lenses Structure and Dynamics (LSD) Survey (Koopmans 
& Treu 2002,2003; Treu & Koopmans 2002a,2003; here- 
after KT02, KT03, TK02a, TK03, or collectively KT). 
The survey takes advantage of the fact that distant early- 
type galaxies are efficient gravitational lenses. By focusing 
on lens galaxies, we can use gravitational lensing analy- 
ses to provide a precise and accurate mass measurement 
(typically 1-5 effective radii R e ), replacing very effectively 



the traditional dynamical tracers at large radii (e.g. X- 
ray; planetary nebulae or globular clusters kinematics) 
that cannot be used in the distant universe. The lens- 
ing analysis is then combined with stellar kinematic mea- 
surements, which provide constraints on the mass distri- 
bution at smaller radii (typically < R c ). The combination 
of the two diagnostics has proved to be very effective (KT) 
since they complement each other: lensing provides a ro- 
bust integrated mass measurement - breaking the mass- 
anisotropy degeneracy of the stellar dynamical analysis - 
after which stellar dynamics provides a handle on the mass 
density profile of the lens. 

The target lenses were selected from the sample of 
known galaxy-scale systems (see e.g. the CASTLES web- 
page at URL http://cfa-www.harvard.edu/castles/) for 
their morphology (E/S0), brightness of the lens (7 < 22) 
and favorable contrast between the lens and the source 
to allow for internal kinematic measurements, and rela- 
tive isolation (e.g. no rich clusters nearby) to simplify as 
much as possible the lens model and reduce the related 
uncertainties. 

Spectroscopic observations using the Keck Telescopes 
have now been completed (a total of nine allocated nights 
between July 2001 and December 2002), yielding exquisite 
internal kinematics for many systems, including 9 early- 
type E/S0 lens galaxies in the range z ~ 0.1-1.0. In 
this paper, we present new data and models for three 
lenses: CFRS.03.1077 (Crampton et al. 2002; Hammer 
et al. 1995; Lilly et al. 1995; ^=0.938, ;z s =2.941 for 
the lens and source respectively), HST1417+5226 (Rat- 
natunga et al. 1995; Crampton et al. 1996; z;=0.810, 
z s =3.399), HST1543+5352 (Ratnatunga, Griffiths & Os- 
trander 1999a; the newly measured redshifts are z;=0.497, 
z s =2.092). We will refer to these lenses as C0302, H1417, 
H1543, respectively. Note that all three objects were 
serendipitously discovered from HST images, i.e. the 
Groth Strip Survey (Groth et al. 1994), the Medium Deep 
Survey (Ratnatunga, Griffith & Ostrander 1999b) and 
the HST Follow-up to the Canada France redshift Sur- 
vey (Brinchmann et al. 1998). Together with the analyses 
of the systems MG2016 and 0047 already presented by 
KT02, TK02 and KT03, this completes the sample of the 
five high- redshift (z <~ 0.5 — 1.0) pressure-supported sys- 
tems targeted by the LSD Survey so far. An analysis of the 
sample properties is presented here. The data and anal- 
ysis of the lower redshift systems and partly rotationally 
supported systems in the current sample will be presented 
in forthcoming papers. 

The paper is organized as follows. In Section 2 we 
present Hubble Space Telescope (HST) imaging and Keck 
spectroscopic observations of the three lens systems. In 
Section 3 we use the photometric and kinematic measure- 
ments to compare the sample of E/S0 lens galaxies to 
the local Fundamental Plane (Djorgovski & Davis 1987; 
Dressier et al. 1987). The offset from the local Funda- 
mental Plane is used to measure the redshift evolution of 
the stellar mass-to-light ratio and thereby constrain their 
star formation history. In Section 4 we present gravita- 
tional lens models of the three lenses, using a modeling 
technique based on a non-parametric source reconstruc- 
tion (e.g. Wallington et al. 1996; Warren & Dye 2003) to 
fully take advantage of the extended nature of the multiply 
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Fig. 1. — Hubble Space Telescope (HST) WFPC2 images in B, V and I band of the three gravitational-lens systems. The right panels show 
a color-composite image. The image sizes are 8"x8" for C0302 and 10"xl0" for H1417 and H1543, and the compass indicates the orientation 
of the images on the sky. Note the bright galaxy G2 (z=0.506) next to the lens in H1543. C0302 is shown off-center because it is close to the 
edge of the WFPC2 chip. 



imaged sources. In Section 5 we introduce two-component 
(luminous plus dark-matter) mass models that will be used 
in Section 6 to perform a joint lensing and dynamics anal- 
ysis and derive limits on the stellar mass-to-light ratio, 
on the inner slope of the dark-matter halo and on the to- 
tal mass density profile of the sample of E/SO galaxies. 
We also consider the complete sample of five high red- 
shift lenses and discuss the evolution of the stellar mass 
to light ratio in terms of stellar population and struc- 
tural/dynamical evolution, the dark-matter mass fraction 
and the limits on the inner slope of the dark-matter halos 
from a joint statistical analysis of the sample. In Section 7 
we discuss the homogeneity of the mass distribution of the 
lens galaxy population, and its implication for lens based 
studies such as the determination of the Hubble constant 
from gravitational time delays, early- type galaxies evolu- 
tion, and the determination of cosmological parameters 
from lens statistics. A final summary is given in Section 8 
and conclusions are drawn in Section 9. 

In the following, we assume that the Hubble constant, 
the matter density, and the cosmological constant are 



H = 65/i65 kms _1 Mpc _1 with h e5 = 1, fl m = 0.3, and 
f^A = 0.7, respectively. Throughout this paper, r is the 
radial coordinate in 3-D space, while R is the radial coor- 
dinate in 2-D projected space. 

2. OBSERVATIONS 

2.1. HST imaging 

Wide Field and Planetary Camera 2 (WFPC2) im- 
ages of the systems CFRS.03.1077, HST1417+5226 and 
HST1543+5352 are available from the HST archive. Five 
exposures each through the F450W filters and the F814W 
filter are available for C0302, with a total exposure time 
of 7000s and 6700s respectively. Four exposures (4400s) 
are available for H1417 through F814W. Sixteen exposures 
are available through filter F606W at 3 different Position 
Angles (PAs), with a total exposure time of 11200s. For 
HST1543, 3 exposures are available through filter F450W 
(8800s), 3 through F606W (9000s) and 2 through F814W 
(6000s). 

The images were reduced using a series of iraf scripts 
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Fig. 2. — Keck-ESI spectra of the E/SO lens galaxies in C0302, H1417 and H1543. A stellar template broadened to the best fitting velocity 
dispersion is overplotted. The bin-sizes along the major axes are indicated in arcseconds. The width of the bins are 1'.'25. 



based on the iraf package drizzle (Fruchter & Hook 
2002), to align the different pointings and perform cosmic 
ray rejection. The images were combined on a O'.'l pixel 
scale. An exception was made for the images of H1417 
through F606W, which were combined in three groups 
according to the HST PA, to avoid problems related to 
distortion correction and complications of the azimuthal 
structure of the Point Spread Function (PSF). The re- 
duced images of the galaxies are shown in Figure 1 to- 
gether with color composite images. Note that for H1417 
only the F606W images at the same PA as for the F814W 
ones are shown and that in these images a bad column of 
WFPC2 runs very close to one of the multiple images, al- 
though this is not the case for the other two sets of F606W 
images (not shown). 

Surface photometry was performed on the F606W and 
F814W images as described in Treu et al. (1999; here- 
after T99) and Treu et al. (2001b; hereafter TOlb). The 
F450W images were not used for surface photometry given 
the low signal to noise of the (red) lens galaxies and the 
large contamination from the (blue) multiple images. The 
galaxy brightness profiles are well represented by an i? 1 / 4 
profile, which we fit - taking the HST point spread func- 
tion (PSF) into account - to obtain the effective radius 
(R e ), the effective surface brightness (SB e ), and the total 
magnitude. The relevant observational quantities of the 
lens galaxies and their errors are listed in Table 2. The 
errors on SB e and R e are tightly correlated and that the 
uncertainty on the combination logi? e — 0.325* B e that en- 
ters the Fundamental Plane (see Section 3) is very small 



(-0.015; see Kelson et al. 2000; TOlb; Bertin, Ciotti & 
del Principe 2002). The restframe photometric quanti- 
ties listed in Table 2 - computed as described in TOlb - 
are corrected for Galactic extinction using E(B-V) from 
Schlegel, Finkbeiner & Davis (1998). 

Astrometry for the system H1417 - the only one where 
the lensed images can be approximated as point images - 
was derived from the two sets of 6 exposures through filter 
F606W that are not affected by a bad column. The two as- 
trometries agree within the uncertainties and are averaged 
to determine the relative offsets between the lens galaxy 
and the multiple images used to constrain the lens model 
in Section 4. 

2.2. Keck Spectroscopy 

The lens galaxies were observed using the Echelle Spec- 
trograph and Imager (ESI; Sheinis et al. 2002) at the 
Keck-II Telescope during four runs on July 21-26 2001, 
February 7-8 2002, June 6 2002, December 7-8 2002. Con- 
ditions were generally photometric with episodes of thin 
cirrus. Between each exposure, we dithered along the slit 
to allow for a better removal of sky residuals in the red 
end of the spectrum. The slit (20" in length) was aligned 
with the major axis of the lens galaxy (C0302, H1417) or 
slightly tilted as to include the massive nearby companion 
(H1543; Figure 1). ESI was used in high resolution mode 
with a ±'.'25 wide slit, yielding a resolution of — 30 km s _1 , 
adequate for measuring the stellar velocity dispersion and 
removing narrow sky emission lines. The centering of the 
lens galaxies in the slit was constantly monitored by means 
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Table 1 
Spectroscopic Observing log. 



Galaxy 


Instrument 


Date 


Seeing 


Exptime 


PA 


C0302 


ESI 


Dec 7,8 2002 


078 


23400s 


110 


H1417 


ESI 


Feb 7 2002 


0'/8 


6300s 


37.2 


H1417 


ESI 


Jul 25,26 2001 


0'/7 


7200s 


37.2 


H1543 


ESI 


Jun 6 2002 


0"6 


14400s 


68 


H1543 


LRIS 


Mar 5 2003 


no 


5400s 


142 



of the ESI viewing camera (all galaxies were bright enough 
to be visible in a few seconds exposure) and we estimate 
the centering perpendicular to the slit to be accurate to 
< O'.'l. Additional details of the observing runs are given 
in Table 1. The ESI data were reduced using the iraf 
package EASI2D 7 as described in KT. The redshifts of the 
lenses are given in Table 2. 

First, the full integrated spectra were used to derive 
the central velocity dispersion of the lens galaxies with 
maximal accuracy. The aperture velocity dispersions are 
measured by comparing broadened stellar templates with 
the observed galaxy spectra in pixel space, as described 
in Treu et al. (1999, 2001) and Koopmans & Treu (2002, 
2003). They are subsequently converted into a central ve- 
locity dispersion (i.e. within a circular aperture of radius 
R e /8), applying an upward correction factor of 1.08±0.04 
(Treu et al. 2001). The central velocity dispersion of G2 - 
the massive companion to H1543 at z=0.506 (see Fig. 1) - 
is found to be o — 263 ± 11 km s^ 1 , by applying the same 
procedure. 



— i — | — i 1 1 — i — | 1 1 — i — i 1 1 — r 
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Wavelength (A) 

Fig. 3.— A Keck-LRIS spectrum (3400-5300 A) of the lensed 
source in H1543. Multiple absorption features clearly identify a 
redshift of z s = 2.092. The lower curve represents the noise level. 

Second, we derive spatially resolved kinematic profiles 
in the following manner. For each galaxy, we define three 
symmetric spectroscopic apertures centered on the lens 

7 developed by D. Sand and T. Treu; Sand et al. (2002,2004) 

8 Using the Gauss-Hermite Pixel Fitting Software, van der Marel 1994. 



galaxy, such that the signal-to-noise ratio of the spec- 
tra integrated within each aperture is sufficient to mea- 
sure 8 a stellar velocity dispersion (see the discussion in 
Koopmans et al. 2003). Portions of the spectra of the 
lens galaxies around the G-band (4304A) for each spec- 
troscopic aperture are shown in Figure 2. The apertures 
are listed in Table 3 together with the measured stellar 
velocity dispersions. Note that our measurement of the 
central velocity dispersion of H1417, a — 224 ± 15 km s _1 , 
is in excellent agreement with the two measurements that 
have been published so far: Ohyama et al. (2002) found 
fj ap = 230 ± 14 km s" 1 within an aperture equivalent to 
a circular aperture of radius 0'.'4; Gebhardt et al. 2003 
measure <7 ap = 202 ± 9 km s _1 within their spectroscopic 
aperture, which they correct to a central velocity disper- 
sion of a = 222 ± 10 km s -1 . 

Since the ESI spectra did not yield the redshift for the 
lensed arc in H1543, additional spectra were obtained us- 
ing LRIS (Oke et al. 1995) on March 5 2003, to exploit 
the favorable contrast at the blue end of the spectrum. 
The l'.'O slit was centered on the arc and aligned with the 
parallactic angle (Table 1). The blue end of the spectrum 
(Figure 3) reveals a set of UV absorption lines typical of 
Lyman-break galaxies (Steidel et al. 1996) that unambigu- 
ously yields the redshift of the arc as z s =2.092 ± 0.001. 

3. THE FUNDAMENTAL PLANE AND THE EVOLUTION OF 
THE STELLAR POPULATIONS 

Early-type galaxies in the local Universe occupy approx- 
imately a plane in the three-dimensional space defined by 
the parameters; effective radius (logi? e ), effective surface 
brightness (SB e ) and central velocity dispersion (log a), 

log R e = a FP log er + /3 FP SB C + 7fp (1) 

known as the Fundamental Plane (hereafter FP; Dressier 
et al. 1987; Djorgovski & Davis 1987). 

Under appropriate assumptions (e.g. Treu et al. 2001), 
the evolution of the FP with redshift can be used to mea- 
sure the star-formation histories of early-type galaxies (e.g. 
Franx 1993). Specifically if we can define an effective mass 
M cx a 2 R e , and if there is no evolution of the slopes «fp 
and /?fp, the evolution of the intercept can be used to mea- 
sure the evolution of the average effective mass to light 
ratio AlogM/L = — 0.4A7fp//?fp (e.g. van Dokkum & 
Franx 1996; Kelson et al. 1997; Bender et al. 1998; Pahre 
1998; van Dokkum et al. 1998; Jorgensen et al. 1999; 
Treu et al. 1999; Ziegler et al. 2001; Treu et al. 2002; 
van Dokkum & Stanford 2003; van Dokkum & Ellis 2003; 

as described, e.g., in Treu et al. 2001. 
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Table 2 

Observed spectro-photometric quantities 



Lens 


C0302 


H1417 


H1543 


redshift (lens) 


0.938±0.001 


0.810±0.001 


0.497±0.001 


redshift (source) 


2.941±0.001 


3.399±0.001 


2.092±0.001 


F814W (mag) 


19.86±0.11 


19.59±0.05 


20.22±0.10 


F606W (mag) 




21.53±0.05 


20.66±0.11 


SB e:F814W (mag/ arcsec 2 ) 


ZZ.o 1 ±U.lo 


Zl. I IzhU.lZ 


on onxn 1 n 


SB e ^ Fe0 ew (mag/arcsec 2 ) 




z4.U0±U.14 


oi on_i_r\ 1 1 
zT.aU±U.ii 


Re,F8i4W (arcsec) 


1.DU±U.1o 


l.Ub±U.Uo 


f\ AA _i_n O/l 
U.41±U.U4 


Re F606W (arcsec) 




1.29±0.13 


0.42±0.04 


b/a=(l-e) 


0.75±0.05 


0.85±0.05 


0.95±0.05 


Major axis P. A. (°) 


-72±5 


34 ±5 


56±5 


(j (kms _i ) 


251±19 


224±15 


116±10 


My - 5 log h 65 (mag) 






-21.40±0.10 


M B - 5 log h 65 (mag) 


-23.30±0.1 


-22.98±0.055 


-20.63±0.11 


R e .y (^-65 lk P c ) 






2.7±0.3 


Re,B (/l^kpc) 


14.1±1.3 


8.6±0.6 


2.8±0.3 


SB e ^y (mag/arcsec 2 ) 






19.34±0.10 


SB C; B (mag/arcsec 2 ) 


20.92±0.14 


20.25±0.12 


20.16±0.11 



Note: The second part of the table lists rest-frame quantities, derived from the observed quantities as described in Section 2. Note that a is 
the central velocity dispersion corrected to a circular aperture of radius R e /8. All quantities in this table assume H = 65 kms" 1 Mpc" 1 , 
fl m = 0.3 and Q A = 0.7. 



Gebhardt et al. 2003; van der Wei ct al. 2004). In the 
following analysis we will adopt as a reference the local 
FP from Bender et al. (1998), i.e. a = 1.25, (5 = 0.32, 
and 7fp = —8.895 — log /150, noting that our results do not 
depend critically on the adopted coefficients of the local 
FP (e.g. Trcu et al. 1999, 2002). 

3.1. Stellar mass-to-light evolution from the FP 

In Fig. 4 we plot the evolution of the effective M/L 
for the five lens galaxies (z w 0.5 — 1.0) analyzed so far 
from the LSD Survey (i.e. 0047-285, MG2016+112 and 
the three systems discussed in Sect. 2), together with the 
published linear fits for the FP M/L evolution of cluster 
and field E/S0 galaxies. The effective M/L evolution for 
the five lens galaxies is d\og(M / L B ) / dz = -0.72 ± 0.10, 
i.e. E/S0 galaxies were on average brighter at z = 1 by 
1.82 ± 0.26 magnitudes in the restframe B band. This 
number assumes that the intercept of the local FP of 
field E/S0 galaxies is the same of that of Coma-Cluster 
galaxies, consistent with the very small environmental 
dependence of the intercept in the local Universe (less 
than 0.1 dex in logM/Ls; e.g. van Dokkum ct al. 
2001; Bernardi et al. 2003). For completeness, we men- 
tion that if the local intercept is allowed to vary, the 
best fit values are dlog(M/L B )/dz = -0.42 ± 0.19 and 
iFP.ienses - lFP,Coma = -0.23 ± 0.14, where the error 
bars are large because of the small number of points and 
the limited redshift range covered by our sample. The 
full LSD sample - including the low redshift objects - is 
needed to simultaneously fit for the local intercept and 
its redshift evolution. In this paper we shall therefore re- 
strict our analysis to the fit of the redshift evolution. We 
note that lens galaxies are not located in rich clusters and 
therefore this sample of lens galaxies is more similar to the 



samples of field E/SOs than those of cluster E/SOs. 

In terms of the passive evolution of a single stellar pop- 
ulation this corresponds to a relatively recent epoch of 
formation z ssp ~ 1.3 (for a Bruzual & Chariot GISSEL96 
population synthesis model with Salpeter IMF and solar 
metallicity; see Treu et al. 2001, 2002 for more details), i.e. 
somewhat younger stars than typically observed for mas- 
sive cluster E/S0 galaxies (for which d\og{M / L b) / dz = 
—0.49 ±0.05 corresponding to z ssp ~ 2; van Dokkum et al. 
1998). However, this evolutionary rate is also consistent 
with a scenario where most of the stars in field E/S0 galax- 
ies are old and formed at z > 2 (for example, MG2016 was 
found to have a star-formation redshift of ^2; KT02) , and 
the relatively fast evolution is driven by secondary bursts 
of starformation contributing <~ 10% to the stellar mass 
between z <~ 1 and today. The latter scenario appears to 
be favored on the grounds of three independent lines of 
evidence: (i) Early-type galaxies are present in the field in 
significant numbers out to well beyond z ~ 1 , inconsistent 
with a sudden creation at z <~ 1.3 (e.g. Treu & Stiavelli 
1999; Chen et al. 2002; Fukugita et al. 2004; Glazebrook et 
al. 2004). (ii) Recent starformation in a fraction of high- 
z E/S0 galaxies is detected from alternative diagnostics, 
such as colors (e.g. Menanteau et al. 2001; van de Ven 
et al. 2003) and absorption and emission lines diagnostics 
(e.g. Trcu et al. 2002; Willis et al. 2002; van Dokkum & 
Ellis 2003). (iii) A "frosting" of younger stars is found 
by comparing detailed stellar population models with the 
spectra of local field E/SOs (e.g. Trager et al. 2000). 

The faster evolution of field versus cluster E/S0 galaxies 
is qualitatively in agreement with the prediction of hier- 
archical models (Kauffmann 1996; Diaferio et al. 2001; 
Benson et al. 2003), however quantitatively the difference 
is smaller than predicted (see Treu 2004 for a recent re- 
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view and discussion of this comparison from a more general 
point of view.) 
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Fig. 4. — Evolution of the FP in rest frame B band. Solid squares 
with error bars are the E/SO lens galaxies analysed in this paper. 
The heavy solid line (KT04) is a linear fit to these points. The other 
lines represent linear fits to measurements from the literature for 
the field FP (Treu ct al. 2002=T02; Gebhardt et al. 2003=G03; van 
Dokkum & Ellis 2003=vDE03; van der Wei et al. 2004=vdW04), the 
cluster FP (vD98=van Dokkum et al. 1998), and the FP based on 
estimates of <tsie using the method of Kochanck ct al. 2000 (Rusin 
et al. 2003=R03; van de Ven et al. 2003=vdV03). Note that the lin- 
ear fits of T02, TK04 (this paper) and vdW04 are almost identical. 
See text for details. 

3.2. Discussion and comparison with previous work 

Our measurement is in good agreement with the re- 
sults - using direct measurements of the central stellar 
velocity dispersion, a - published by Treu et al. (2002) 
who find d\og(M/L B )/dz = -0.71±q\1 and by van der 
Wei et al. (2004) who find -0.71 ± 0.20. In contrast, 
van Dokkum & Ellis (2003) find a marginally slower 
evolution dlog(M / L B ) / dz = -0.55 ± 0.05, while Geb- 
hardt et al. (2003) measure a much faster evolution, i.e. 
2.4 magnitudes to z = 1, corresponding to a value of 
d\og{M / L b) I dz w —0.96. At face value these results 
appear inconsistent at the 1-2 a level. Assuming that 
mutually consistent measurement techniques have been 
adopted 9 , the differences could arise for a variety of rea- 
sons. 

A first possible explanation is that the various samples 
cannot be directly compared because of subtle differences 
between the morphological classification schemes. Indeed 
van Dokkum & Ellis (2003) found that the 2 overluminous 
E/S0 galaxies in their sample of 9 showed asymmetric fea- 
tures in the Hubble Deep Field images. A larger sample 
of objects with deep imaging is necessary to quantify this 
effect. 



Second, luminosity-selected samples favor overluminous 
objects and are therefore biased toward faster evolution. 
The extent of the bias depends on the intrinsic scatter of 
the FP and on the magnitude limit. It can be corrected by 
taking into account the selection procedure in the analysis 
(Treu et al. 2001, 2002). To further complicate matters, 
the evolutionary rate could be a function of mass, with 
more massive galaxies evolving slower than less massive 
ones (e.g. van der Wei et al. 2004; note also that H1543, 
the lens galaxy with the largest offset from the local FP is 
also the least massive one), resulting in a change of the FP 
slopes with redshift. In this case, the mean evolutionary 
rate would be a function of the mass range of the sample 
(Treu et al. 2002). Unfortunately, luminosity selection 
could also mimic a change in the slope, because less mas- 
sive galaxies would make the cut only if they are overlumi- 
nous (Kelson et al. 2000). We note here that a lens sample, 
although in principle mass selected, could suffer from the 
same kind of bias because of the effective magnitude limit 
imposed by spectroscopic follow-up. The statistical anal- 
ysis of larger samples of lens galaxies is essential to deter- 
mine simultaneously the evolution of the intercept, slopes 
and scatter of the FP, correcting for selection effects. 

A third possible explanation is cosmic variance. If in- 
deed the secondary burst scenario is correct, at any given 
time the majority of E/S0 galaxies would be observed to 
follow a quiescent passive evolution path in the FP space, 
while a fraction of E/S0 galaxies would be observed within 
1-2 Gyrs after the secondary burst, while overluminous. 
The fraction of overluminous E/S0 galaxies would depend 
on the duty-cycle of secondary bursts. If, for example, 
~10% of the stellar mass is formed in each secondary burst, 
each galaxy undergoes one secondary burst between z = 1 
and z = 0, and the bursts are detectable for 2 Gyrs, then 
a quarter of the E/S0 galaxies between < z < 1 would 
be overluminous. Therefore, only a handful of overlumi- 
nous galaxies would be observed in the current samples of 
galaxies and small number statistics could be dominating 
the uncertainties. 

Following Kochanck ct al. (2000), both Rusin ct al. 
(2003) and van de Ven ct al. (2003) recently used image 
separations (or Einstein Radii) of arcsecond-scale strong 
lens systems to estimate the central velocity dispersion of 
E/S0 lens galaxies and construct a FP. The key assump- 
tion is that lens galaxies have isothermal mass profiles 
(i.e. p oc r~ 2 ) and therefore the central velocity disper- 
sion can be obtained directly from the image separation 
with the assumption that a w ctsie (see Kochanek et al. 
2000). Under this assumption, Rusin et al. (2003) find 
d\og{M/L B )/dz = -0.56 ± 0.04, and van de Ven et al. 
(2003) - using the photometry and image separations from 
Rusin ct al. (2003) - find d\og(M / L B ) / dz = -0.62±0.13, 
using larger errors and a different weighting scheme. The 
fact that these estimates are so close to the direct mea- 
surements discussed above is indeed remarkable. Not only 
do many of the above arguments, related to selection ef- 
fects and small sample statistics, also apply to lens galaxy 
samples, departures from isothermal mass profiles, or ef- 
fects such as a mass-sheet degeneracy, introduce additional 
sources of uncertainty. In other words, the relatively good 
agreement between the direct and indirect methods tells 



As suggested by the good agreement between the FP parameters measured by different groups, see Section 2.2 and Gebhardt et al. 2003 
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Table 3 

Kinematic data along the major axis of the lenses 
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Note: The adjacent rectangular apertures are indicated, as well as the measured aperture velocity dispersions (ex), their uncertainty (Act), and 
the average S/N per A in the region used for the kinematic fit. 



us that the isothermal approximation is probably not dra- 
matically wrong. 

In the next sections, via a joint lensing and dynam- 
ics analysis, we will further examine the accuracy of this 
approximation. However, before we proceed to the full 
analysis we can gain some insight by looking at the three 
objects in our sample that are also in R03 and vdV03: 
0047, H1417, MG2016. The offset from the local FP 
are in agreement within the errors for 0047 and MG2016 
(see also R03), while for H1417 we measure <~ 0.1 dex 
more evolution for log(M/LB). The difference is entirely 
due to the difference between our direct measurement of 
the stellar velocity dispersion (a = 224 ± 15km s _1 ) and 
the velocity dispersion inferred from the image separation 
(~290km s _1 ; vdV03, R03; next Section). Hence, in at 
least one case, the image separation underestimates the 
stellar velocity dispersion. As we will show in Section 7, 
this is the case in three out of five lens systems in our 
sample. 

4. GRAVITATIONAL LENS MODELS 

As in previous work on LSD lens systems (KT), we 
model the three systems (C0302, H1417 and H1543) with 
a Singular Isothermal Ellipsoid (SIE) mass distribution: 



Ksm(x,y) 



q 2 x' 2 



(2) 



2vY^ 

with bi = 47t(<tsie/c) 2 (Dd Dd s /D s ) and the major axis 
aligned north-south (Kormann et al. 1994). We define 
(x',y') (in radians) in Eq.2 as the frame centered on the 
lens (xi,yi) and aligned with the PA of the lens 9i. Be- 
cause the mass enclosed by the elliptical critical curve 
is independent of the flattening of the mass distribution 
q = {b/a), we find that i?Einst = D^bi corresponds to 
the equivalent Singular Isothermal Sphere (SIS) Einstein 
radius. The equivalent SIS mass (Msinst) is that en- 
closed by the critical curve. Both i?Einst and Msinst are 
needed in the joint lensing and dynamical analysis. We 
refer to KT03 for additional discussion of the models. We 



note that the enclosed mass Msinst is nearly independent 
of the choice of mass model and depends predominantly 
on the monopole of the lens potential (Kochanek 1991). 
We also allow for an external shear with strength 7 ext 
and position angle # C xt, corresponding to the potential 
4>ext(x, y) = -R 2 (7ext/2) cos 2(6 - ext ), where we define 
R 2 = (x — xi) 2 + (y — yi) 2 to be the square of the distance 
to the lens center (see Keeton et al. 2002). We emphasize 
that our lens model is fully two-dimensional and there- 
fore takes correctly into account deviations from circular 
symmetry (e.g. Sand et al. 2003, Dalai & Keeton 2003, 
Bartelmann & Meneghetti 2003). 

Only H1417 has well-defined lensed images and fluxes 
(Fig. 1). For this lens system, previously modeled by 
Knudson, Ratnatunga & Griffiths (2001), we use the 'tra- 
ditional' modeling technique where the source is repre- 
sented by a single point in the source plane with a given 
flux, which is magnified and mapped onto the lensed im- 
ages in the image plane. We use the code described in 
Koopmans et al. (2003). To model the extended images of 
C0302 and H1543, we implemented a code (see Appendix) 
that incorporates techniques described in Wallington et al. 
(1996) and Warren & Dye (2003), allowing the full use of 
all the lensing information contained in the images. 

4.1. HST14176+5226 

We first model the lens system with a SIE mass distri- 
bution and external shear only. The lens strength (bi), 
ccntroid (xi,yi), ellipticity (q) and position angle (9i) are 
free parameters, as well as the external shear, the source 
position and its flux. The lens-ccntroid and position angle 
are further constrained by the observational priors given 
in Tables 2 and 4.1. The best model has x 2 /NDF = 38 
for NDF=4, using the image constraints in Table 4.1, and 
is clearly not satisfactory 10 . 

The next order of observable complexity that can be 
introduced is a linear gradient in the surface density 
of the mass model and has a potential of ip g (x,y) = 
(i? 3 / 4 ) |Vk s | cos(0- 9 g ) (see Keeton 2002). This gradient 



10 The x 2 value is dominated by the lens position (x! 2 =84) and the lensed image positions (xf=65). The flux-ratios contribute only Xr = 2- 
Since the images are somewhat extended (see Fig. f ) and have no apparent compact structure that could be affected by stars (i.e. microlcnsing) 
or substructure, one would not expect significant flux-ratio anomalies. 
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can be the result of an external perturber (e.g. group or 
cluster), but also an internal asymmetry of the lens (e.g. 
M=l mode). We assume n g = for a line through the 
lens centroid, i.e. R 2 = (x — xi) 2 + (y — yi) 2 , such that 
the convergence gradient adds no mass inside a symmetric 
aperture on the lens centroid (xi,yi). This adds two more 
free parameters (i.e. |V« g | and 6 g ). The best model in this 
case has % 2 /NDF = 3.2 for NDF=2, which is considerably 
better than a model without a convergence gradient. We 
feel that adding more free parameters is no longer justified, 
since both % 2 and NDF arc small. 



Object 


ARA 


A DEC 


r 


A 


-1.266±0.006 


-1.139±0.006 


0.81±0.16 


B 


-0.843±0.006 


+0.918±0.006 


0.65±0.08 


C 


+0.792±0.006 


+1.321±0.006 


1.00±0.13 


D 


+0.814±0.006 


-0.803±0.006 


0.57±0.18 


G 


+0.000±0.004 


+0.000±0.004 





Note: The column r lists the F814W flux ratios of the multiple 
images normalized to image C, taken from the CASTLcS web 
database. 



We notice that the gradient points in the direction of 
the major axis of the SIE mass distribution and not in 
the external shear direction. The direction of the gradi- 
ent and the major axis of the SIE coincide within 11° and 
are consistent at the 2.5 a level. We also notice that the 
agreement improves to within ^9° if the major axis is not 
constrained by a prior. To further examine this possible 
alignment, we also tested the SIE plus shear and gradient 
model on PG1115+080 (a model similar to TK02b). We 
find that for the best model, the position angle of the gra- 
dient and external shear agree, both pointing to the con- 
firmed external perturber, which is a compact group ~13" 
from the lens (Kundic et al. 1997). In the case of H1417, 
however, there exists no obvious group or cluster in the 
~ 1' field around the lens system 11 and the gradient and 
shear position angle differ at the 12 a level, nor are any of 
the galaxies around the lens massive enough to account for 
the observed gradient of |V« g | = 0.102 ± 0.015. A simple 
argument based on a SIS perturber shows that one would 
expect its distance from H1417 to be related to the Ein- 
stein radius of the perturber by (i?/i?Einst) 2 ~ l/(2|V/c g |), 
or R s=a 2.2i?Einst- A massive cluster with i?Einst = 30" 
would still be detectable in the field around the lens. A 
single massive galaxy with i?Einst = 2" would be within 
R w 5" from the lens - in the direction of increasing « g , 
i.e. 8 g w 200°. Neither is obviously found. We therefore 
conclude that the detected gradient is most likely associ- 
ated with an asymmetry in the lens mass distribution. A 
significant gradient is also required for models with density 
slopes other than isothermal. 

The best SIE model is shown in Fig. 5 and all model 
parameters are listed in Table 4.1. The equivalent SIS ve- 
locity dispersion, mass and Einstein radius of the lens are, 
ctsie = 290±8km s" 1 , M Ei „ s t = (70.8±7.6) • 10 10 M and 
ifeinst = H.4 ± 0.6 kpc (l?41±(y/08), respectively (68% 
CL errors). Our derived Einstein radius (l'/41) agrees to 
within <0'/01 with that derived by Rusin et al. (2003). As 
noted the gradient adds no mass. 

To further test the robustness of the enclosed mass mea- 
surement, we examine three more sources of uncertainty. 
First, we find that the mass enclosed by a circular aperture 

11 Although an overdensity of galaxies at z ~ 0.8 in the Groth Strip 



with l'/41 radius, gives a mass that is only 2.0% lower than 
our best estimate of the mass within the elliptical critical 
curve. Second, the best % 2 models with density slopes 
that arc 40% steeper/shallower than isothermal have en- 
closed masses different by only — 2.0%/+1.3%. Third, if n g 
is caused by an external perturber, a mass-sheet (K s heet) 
must be associated with it, since truly negative values of 
k„ arc not allowed. An estimate of the external convcr- 
gence can be obtained from the external shear, assuming 
an isothermal mass distribution, i.e. K s hcct ~ Text- If we 
set Ksheet = 7ext 5 we find that the velocity dispersion, en- 
closed mass and Einstein radius decrease by 2.5%, 10% and 
5%, respectively. However, in the absence of any evidence 
for external perturbers that could result in a non-local 
convergence, we conclude that the mass measurement is 
robust, reliable and does not introduce a bias in the final 
lensing plus dynamics analysis (Section 6). To account for 
most of the uncertainties we assume a total 1-a errors of 
10% on the enclosed mass. 

4.2. CFRS0S.1077 

The lensed arc and counter image in the system C0302 
(Crampton et al. 2002; Hammer et al. 1995; Lilly et al. 
1995) do not have enough structure to allow a simple one- 
to-one mapping between them. To determine the total 
enclosed mass, we implemented a lens code that combines 
elements from Wallington et al. (1996; W96 from hereon) 
and Warren & Dye (2003; WD03) and can model extended 
lensed images on a grid (i.e. CCD image). In the Appendix 
we outline the general features of the code and where it 
differs from W96 and WD03. 

In the discovery paper, Crampton et al. (2002) mod- 
eled the lens using an NIE mass model with non-zero core 
radius. They find a good fit to the arc and counter im- 
age, using an axisymmetric source model just outside the 
cusp. A high velocity dispersion of 387±5 km s _1 is in- 
ferred from their lens model, leading to the conclusion that 
the lens galaxy is as faint as present-day ellipticals of sim- 
ilar central velocity dispersion and thus shows no sign of 
passive evolution. Since there are several dwarf-like galax- 
ies near the main lens galaxy (Fig. 1), possibly indicating 
the presence of a small group, the apparent underlumi- 
nous nature of the lens galaxy might be a result of an 
increased dark-matter fraction inside the Einstein radius 
from a group halo. Groups have been found near a number 
of other lens systems, hence they might not be uncommon 
and should be accounted for in lens models where neces- 
sary (e.g. Lehar et al. 1997; Kundic et al. 1997a&b; Tonry 
1998; Tonry & Kochanck 1999, 2000; Blandford, Kundic 
& Surpi 2001; Keeton, Christlcin & Zabludoff 2000; Fass- 
nacht & Lubin 2002; Johnston et al. 2003). 

In our model, we associate the mass centroid with the 
light centroid of the lens galaxy and model its mass dis- 
tribution as an SIE with the lens strength, position angle 
and ellipticity as free parameters. The position angle is left 
free (as opposed to H1417) because the nearby companions 
(Fig. 1) could introduce a difference between the luminous 
and mass position angle. External shear is added. 

The source grid is 39x39 pixels (40 mas per pixel). Pix- 
els in the image plane within 1-a of the noise level are 
clipped and do not participate in the determination of 

t be present; Koo et al. 1996; Im et al. 2002. 
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Fig. 5. — Best-fitting SIE lens model of H1417. The left panel shows the position of the lens (solid circle) and the multiple images (solid 
squares) on the image plane, together with the critical line (thick solid line). The caustics (thin lines) and source position (solid star) on the 
source plane are also shown. The right panel shows the time delay surface with constant time-delay contours increasing from 0.0 (at A) to 83 
days in steps of 8.3 days. 




Fig. 6. — The best LENSGRID (Section 4) reconstructions of the extended arcs and counter images in C0302 (left) and H1543 (right). Panels 
(a) show the original image with lens-galaxy subtracted. Panels (b) show the best reconstruction. Panels (c) show the difference between 
observation and model and panels (d) show the source in the source plane, regridded and smoothed by a Gaussian with a FWHM=0'.'08 to 
highlight its structure. The curves are the critical (upper-left) and caustic (lower-left) curves, respectively. 
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Fig. 7. — Likelihood contours for a joint lcnsing and dynamical analysis of C0302 (upper row), H1417 (middle) and H1543 (lower). Panel 
(a): The contours in the M+fLg-'y plane for an isotropic mass model, with (thick) and without (dashed) the FP constraints. Panel (b): as 
(a) for a radially anisotropic model. Panel (c): contours in the ri vs 7' plane. Panel (d): contours in the 7' — &i so plane. Note that the 
effective slope 7' changes very little for bi so = — > — 1. 



the best model (clipping at 2-3 a results in a very similar 
model). Image-plane pixels that map outside the source 
grid and visa versa arc masked. We minimize the re- 
duced x 2 , to properly account for the variable number 
of participating source and image-grid pixels. We use an 
appropriate PSF generated with Tiny-Tim. A value of 
A=0.02 for the smoothing parameter (see Appendix) leads 
to x 2 /NDF w 1. 

Since the external shear is large, we find a strong degen- 
eracy between the lens-galaxy ellipticity and position an- 
gle and the external shear. The constraint on q = (b/a) is 
therefore relatively weak. Since the external shear appears 



to align with the position angle of the light distribution, 
not the mass model, we suspect that some of the external 
shear is in fact 'internal' shear due the stellar component 
and that the lens galaxy is embedded in a larger misaligned 
structure (e.g. a group halo). The presence of at least five 
small dwarfs around C0320 would support this (Fig. 1). 
These degeneracies, however, have negligible effect on the 
determination of enclosed mass within the images. 

The best SIE plus shear model is shown in the left panels 
of Fig. 6 and the mass- model parameters are listed in Ta- 
ble 4.3. Wc have refrained from calculating precise formal 
errors - as for H1417 - which is extremely difficult given 
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the variability of NDF (i.e. the number of participating 
pixels) and the free choice of source-pixel size and A. It 
is also computationally very expensive if the mass mod- 
els is not fixed (see WD03 for the case when it is fixed). 
We plan to further refine the code - currently written in 
IDL -, increase its speed and allow for a full non-linear 
error analysis. Nevertheless, a conservative upper limit 
of <5% can be set on the 1-cr error on the Einstein ra- 
dius, being roughly the width of the arc divided by its 
distance to the lens centroid. The equivalent SIS velocity 
dispersion, mass and Einstein radius of the lens in C0302 
from the best model then become, <tsie = 294 ±8 km s -1 , 
M Einst = (67.0±6.7)-10 10 M and i? E i„st = 10.6 ±0.5 kpc 
(l'/24±0'/06), respectively. 

The dominant component of the source straddles the 
cusp on the ouside and appears compact and relatively 
symmetric with possibly some indication of extended 
structure around it. The shape and position in the source 
plane are similar to those found by Crampton et al. (2002), 
even though their stellar velocity dispersion, from a non- 
singular model, is much higher than our SIE velocity dis- 
persion. This appears to be a result of the use of a non-zero 
core radius and their definition of a, which increases with 
increasing core radius (Francois Hammer, private commu- 
nication). We note that our definition of ctsie i s m accor- 
dance with previous work, including that of Kochanek et 
al. (2000) and Rusin et al. (2003). 

4.3. HST15433+5352 

The modeling of H1543 proceeds very similar to that of 
C0302, with the following differences. First and foremost, 
there are several massive nearby perturbers that have to 
be accounted for in the model. The strongest is a mas- 
sive galaxy (G2) approximately 4'.'7 to the east of the lens 
galaxy (Gl) with a measured central stellar velocity dis- 
persion of & — 263 ± 11 km s _1 . We model this galaxy as 
a SIS (Eq.2 with q = 1.0) with ctsie = &■ 

Another galaxy (G3) at the same redshift of Gl and G2 
was serendipitously detected in the LRIS slit, 18" away 
from Gl on the opposite side of G2. Since H1543 falls 
near the edge of the WFPC2 field, no HST images are 
available of G3. However, images of the region are avail- 
able from the Sloan Digital Sky Survey (SDSS) Data Re- 
lease 1 (Abazajian et al. 2003). A visual inspection of the 
images shows indeed a galaxy at the location of G3, and 
other galaxies with similar colors and luminosities in the 
vicinity, consistent with the presence of a group (see the 
discussion on groups for C0302). 

Our best SIE plus shear model is shown in Fig. 6. Rel- 
atively little regularization is needed (A ~ 0.003). The 
source is relatively compact, which lends further credit to 
the model. As expected from the presence of the nearby 
aligned perturbers (G2 and group), we find the external 
shear, -f cxt = 0.17, to be large and dominant over the mass 
ellipticity of lens galaxy. To avoid degeneracies between 
external shear and lens galaxy ellipticity, we therefore re- 
stricted ourselves to a SIE mass model with q = 0.95 fixed. 
We also tested SIE plus shear models with varying ellip- 
ticities and find no strong differences between their criti- 
cal and caustic structures, nor between the enclosed mass 
within the equivalent SIS Einstein radius. The shear po- 
sition angle (^72°) aligns nearly perfectly with the line 



between Gl, G2 and the compact group (^67°), suggest- 
ing the external shear is real and most likely due to the 
group. The shear of G2 is already accounted for by its SIS 
mass model. 

Table 4 
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Note: The lens centers of C0302 and H1543 are fixed at the 
observed galaxy centroids. The sky PA values are given. 

The equivalent SIS velocity dispersion, mass and Ein- 
stein radius of the lens in H1543 from the best model are, 
asm = 139 ± 7 km s" 1 , Af E inst = (3.4 ± 0.7) • 10 10 M 
and i?Einst = 2.4 ± 0.4 kpc (0'.'36±0'.'04), respectively. As 
for H1417 we can estimate the mass associated with the 
external mass distribution from the external shear. If a 
mass sheet with K s h cc t ~ 7oxt ~ 0.17 contributed to the 
image separation, one would overestimate the mass of the 
galaxy by ^17%. To account for this uncertainty, we adopt 
a conservative error of 20% on MEinst , twice that of the 
other two systems. This range also covers the majority of 
models, when varying different assumptions in the models 
(e.g. ellipticity). 

We note that if G2 and the group are dynamically asso- 
ciated with the lens galaxy, their dark-matter mass halos 
also contribute to the mass inside the Einstein radius of 
the lens (Gl), affecting both lensing and stellar dynamics. 
Therefore, one should not regard this as a systematic effect 
that should be removed like a mass-sheet, since this mass 
truly contributes to the inner slope of the dark-matter halo 
of the lens galaxy. We discuss this important point in more 
detail in Sect. 7. 

5. THE TWO-COMPONENT MASS MODEL 

Following our previous papers (TK02, KT03), in the 
lensing plus dynamics analysis, we model the lens galaxies 
as a superposition of two spherical components, one for the 
luminous stellar matter and one for the dark-matter halo. 
The luminous mass distribution is described by either a 
Hernquist (1990) 

< \ M * r * ei\ 

Plum(r) = -— Tg (3) 

or a Jaffe (1983) model where M* is the total stellar mass. 
For consistency with TK02 we will show primarily the 
results obtained with the Hernquist profile and discuss 
how they change using a Jaffe profile where relevant. The 
dark-matter halo is modeled as: 



Pbm,o r h 



(4) 



r l ( r 2 +r 2)(3- 7 )/2 

which closely resembles an NFW profile for 7 = 1 , and has 
the typical asymptotic behavior at large radii found from 
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numerical simulations of dark-matter halos oc r~ 3 (e.g. 
Ghigna et al. 2000). In accordance with the CDM picture 
(e.g. Bullock et al. 2001) we expect the break radius 
to be much larger than the effective and Einstein radii. 
Therefore, in the following we will set ^> i?Einst, effec- 
tively equivalent to oo, i.e. we describe the dark- matter 
halo as a simple power-law pdm oc r~ 7 in the region of 
interest. To further explore the effects of the distribution 
of mass at large radii we have done tests with (i) a dark- 



matter halo, as in Eq. 4, but falling off as 



at large 



radii (equivalent to a Hernquist or Jaffe model for 7 = 1,2 
respectively) and (ii) values of r& as small as < R c . In all 
cases we find that the effects on the stellar velocity dis- 
persion due to changes in the break radius and outer slope 
are negligible if > 3i? . In the current ACDM mod- 
els break radii as small as r& < 3i? are highly unlikely 
for most galaxies, since R c is typically only several kpc. 
Our approximation of r& — > 00 is therefore justified and 
the resulting constraints on 7 can be compared with the 
inner dark- matter mass slope from simulations after bary- 
onic collapse (e.g. cooling), which can alter the inner slope 
of the dark-matter halo but much less so the outer slope. 

We assume an Osipkov-Merritt (Osipkov 1979; Merritt 
1985a, b) parametrization of the anisotropy of the stellar 
mass distribution or a constant (3(r) model 



1 - 



n > 



(5) 



6 iso e [-i,+i] 



where ag and ay are the tangential and radial components 
of the velocity dispersion and r, is called the anisotropy 
radius. Note that j3 > by definition for the OM model, 
not allowing for tangcntially anisotropic models. At infi- 
nite radii, Osipkov-Merrit models become completely ra- 
dial. Although this behavior is not commonly found within 
the inner regions of E/S0 galaxies probed by observations 
(e. g. Gerhard et al. 2001; see also van Albada 1982 and 
Bertin & Stiavelli 1993 for theoretical grounds), it has lit- 
tle impact in the case considered here, since the pressure 
tensor only becomes significantly radial well outside the 
Einstein Radius and in projection is significantly down- 
weighted by the rapidly falling luminosity-density profile. 
To test tangential anisotropy, we also considered models 
with constant anisotropy [3{r) = fe; so varying from —1 to 
+1. Whereas for b lso = 1 — > the behavior is very similar 
to ri = — > 00, the effect on the infered mass slope 7' for 
(3 = — 1 to is almost negligible (panels (d) in Fig. 7). We 
will therefore not further consider these models in detail, 
but only mention them when necessary. 
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Fig. 8. — Upper panel: The effective slope as a function of the ra- 
tio between the Einstein Radius and the effective radius; solid points 
are obtained for isotropic models (r» = 00), open points for radially 
anisotropic models (ri = R e )- Lower panel: The ratio between cen- 
tral velocity dispersion (a) and velocity dispersion of the best fitting 
SIE model as a function of the ratio between the Einstein Radius 
and the effective radius. The horizontal dashed lines represent the 
expected values for isothermal lens models. The measured values 
(isotropic models only) for PG1115 (TK02b) and B1608 (K03) - 
not included in the LSD sample - are also shown for completeness 
as open triangles. See Sections 6 and 7 for discussion. 

The line-of-sight velocity dispersion is obtained by solv- 
ing the three dimensional spherical Jeans equation for the 
luminous component in the total gravitational potential 
and computing the luminosity-weighted average along the 
line of sight (e.g. Binney & Tremaine 1987; Kochanek 
1994). We correct for the average seeing during the ob- 
servations, and average the velocity dispersion - weighted 
by the surface brightness - inside the appropriate rectan- 
gular apertures. For completeness, we rescale the model 
apertures such that their projection on the axisymmet- 
ric model is equivalent to their projection on an elliptical 
galaxy with an axial ratio of b/a. This has minimal ef- 
fects on the model velocity dispersions - much smaller (i.e. 
< 1%) than the observational errors - since the observed 
and model dispersion profiles are typically very flat. The 
uncertainties on seeing, aperture size, and galaxy centering 
are taken into account as systematic errors in the following 
discussion. Additional discussion of our mass profiles and 
dynamical model can be found in TK02 and KT03. 

6. A JOINT LENSING AND DYNAMICAL ANALYSIS 

We are now in the position to use the measurements 
derived in Sections 2.1, 2.2 and 4 to constrain the free pa- 
rameters in our two component mass models. Surface pho- 
tometry gives directly r* = i? c /1.8153 for the Hernquist 
model, assuming that the stellar mass-to-light ratio is con- 
stant. The mass enclosed by the Einstein radius is used 
to obtain Pdm,o, given the other parameters. Likelihood 
contours of the three remaining parameters (i.e. M*/Lb, 
7 and r{) are then obtained by comparing the velocity dis- 
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pcrsion profiles from the models with the observed ones. 

6.1. The FP as an additional constraint 

An additional constraint can be obtained using the off- 
set of the galaxy from the local FP by introducing one 
further assumption. If the evolution of the effective mass 
to light ratio A log(M/L B ) is equal to the evolution of the 
stellar mass-to-light ratio Alog(M*/LB), the stellar mass- 
to-light at redshift z is related to the stellar mass-to-light 
ratio at z = by 




where the first term on the right hand side of the equa- 
tion can be measured for local E/SO galaxies. Using the 
local value of (M*/L B )o = (7.3±2.1) h 65 M /L Bi0 , deter- 
mined from data by Gerhard et al. (2001) as in TK02, we 
infer M*/L B = (1.9±0.5) h 65 M /L s , M*/L B = (1.6± 
0.4) h 65 M /L B , , M*/L B = (2.1 ±0.6)^66 M /L Bj0 for 
the lens galaxies in C0302, H1417 and H1543, respectively. 
In the next sections we will compare these measurements 
with the independent ones obtained from the joint lensing 
and dynamics analysis, finding a good agreement. Then, 
we will also use the FP measurements as a further con- 
strain to the two-component models. However - since this 
determination relies on a non-trivial assumption (Eq. 6) - 
we will present both the results that include this constraint 
and those that do not. 

6.2. Single power-law mass models 

Before considering the full two-component models (Sec- 
tion. 5) let us first consider a simplified family of models 
to explore the properties of the total mass distribution, 
since this is of particular relevance to studies of, for ex- 
ample, the value of H from time-delays and also lensing 
statistics. 

As in TK02 and KT03 this family of models consists of 
a total luminous plus dark-matter mass distribution that 
follows a single power-law p t ot °t r 7 within the region 
of interest, where 7' is called the effective slope. Hence, 
the luminous mass is assumed to be a trace component 
in the potential, with M*/Lb = 0. The two remaining 
free parameters, 7' and ?*j are constrained with the veloc- 
ity dispersion profile, yielding the results shown in panels 
(c) of Fig. 7. In panel (d), the constant (3 models are 
shown, displaying a similar behavior for 6i so > 0, whereas 
for &i so < (tangential anisotropy) almost no effect is seen 
on the value of 7'. 

The best-fit values of 7' depend on the anisotropy of 
the velocity ellipsoid. As expected, an isotropic velocity 
ellipsoid (rj = 00 or 6 iso = 0) leads to a larger value of 7', 
whereas for the more radial orbital structures (r, — > 0), 
a smaller value of 7' is needed. The most likely values of 
7' for two representative cases (r» = 00 and r» = R e ) are 
listed in Table 5. The corresponding values for MG2016 
and 0047 (TK02a, KT03) are also listed for completeness. 

The results are also shown in Figure 8, where we plot 
7' as a function of i?Einst/-Rc- The average slope from the 
five systems in our sample is (7') = 1.75±0.09 with a large 
rms of 0.20 (isotropic), or (7') = 1.57 ± 0.16 with an rms 
of 0.35 (anisotropic). Extremely radial orbits (7^ < R e ) 



can probably be ruled out, both on observational grounds 
(e.g. Gerhard et al. 2001) and on theoretical grounds, since 
they would lead to instabilities (Merritt & Aguilar 1985; 
Stiavelli & Sparke 1991), whereas tangential anisotropy 
can not be ruled out, it has a negligible effect (see Panels 
(d) in Figure 7). 

In the lower panel of Figure 8, we also show /sie, i-e. 
the ratio between the central velocity dispersion and the 
velocity dispersion of the best fitting Singular Isothermal 
Ellipsoid. This number is independent of the choice for 
the dynamical model. The average is (/sie) = 0.87 ± 0.04 
with an rms of 0.08, lower on average than that based on 
the expectation that a w ctsie (Kochanek 1994; Kochanek 
et al. 2000). 

6.3. Luminous and dark-matter mass decomposition 

Let us now consider the two component mass models. 
Once again, we examine the two cases of r, = 00 and 
Ti = R c , delineating a conservative range of physical mod- 
els, and derive likelihood contours in the M^/L B ~J plane. 
The likelihood contours are shown in the panels (a) and 
(b) of Fig. 7. The dashed lines represent the contours ob- 
tained without the FP constraint on M*/L B , while the 
solid contours include the FP constraint (Sect. 6). The 
main effect of including the FP constraint is to rule out 
regions with low stellar mass-to- light ratios. In general, 
the shape of the contours is well understood: the outer 
luminous component is on average steeper than the dark- 
matter component and therefore smaller values of M*/L B 
require larger values of 7, to compensate and produce a to- 
tal mass profile as steep as required by the kinematic data. 
Increasing the radial anisotropy implies a more shallow to- 
tal mass profile and therefore a smaller value of 7. 

In this subsection we discuss constraints on the amount 
of dark matter within the Einstein radius, through mea- 
surements of the total and stellar mass to light ratio. In 
the next subsection (Sect. 6.4), we will use these values to 
determine the cosmic evolution of the stellar mass to light 
ratio and hence the star formation history of E/S0 galax- 
ies. In light of these two goals we will determine the frac- 
tion of dark matter marginalizing the likelihood contours 
shown in Figure 7 over 7 in two ways: (i) Including the 
FP constraint as a prior, to obtain the most precise mea- 
surement of M*/L B . (ii) Not including the FP constraint, 
to obtain an independent measurement of the cosmic evo- 
lution of M*/L B . The results are listed in Table 6. The 
68% confidence limits - around the maximum-likelihood 
value - are uniquely determined from the probability dis- 
tribution function of M^/Lb (after marginalising over 7) 
by the two values of M*/L B that have equal probability 
densities and enclose 0.68 of the probability. Since for low 
values of M*/L B , the probability density as function of 7 
is only large for a small range around 7 ^ 2, it is nearly 
constant over a large range of 7 for M*/ L B ~ 2—3 (Fig. 7). 
Hence, after marginalisation, M*/L B has a clearly-defined 
maximum- likelihood value and 68% confidence limits. 

The total mass to light ratio enclosed within the Ein- 
stein radius (M tot /L B )(< i?Einst) is also listed in Table 6 
for comparison. Note that (M to t/ L B )(< ^Einst) is consid- 
erably higher than the limit on M^jL B for all five E/S0 
galaxies, implying that the galaxies cannot be described by 
a constant mass-to-light ratio model. Hence, the presence 
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Table 5 

Summary of lens/dynamical model results - I 



Galaxy 


7' (iso) 


7' (aniso) 


CTSIE 

(km s" 1 ) 


/sie 


R-Einst 

(arcsec) 


Re 

(arcsec) 


0047 

C0302 

H1417 

H1543 

MG2016 


1.90±0.05±0.1 
1.70±0.05±0.1 
1.75±0.05±0.1 
1.40±0.20±0.2 
2.00±0.10±0.1 


1.7±0.05±0.10 
1.3±0.10±0.10 
1.7±0.05±0.10 
1.15±0.05±0.10 
2.0±0.10±0.10 


252±4 
294±8 
290±8 
139±7 
331±10 


0.91±0.06 
0.85±0.06 
0.77±0.05 
0.83±0.13 
0.99±0.10 


1.34±0.01 
1.24±0.06 
1.41±0.08 
0.36±0.04 
1.56±0.02 


0.82±0.12 
1.60±0.15 
1.06±0.08 
0.41±0.04 
0.31±0.06 



of a mass component spatially more extended than the 
luminous component is required. We identify this compo- 
nent with the dark-matter halo. Quantitatively, we find 
that for all five E/SO galaxies (including MG2016 and 
0047) , the no-dark-matter-halo scenario is excluded at the 
> 99% C.L. In other words, a constant mass-to-light ra- 
tio model is too steep to satisfy simultaneously the lensing 
and dynamical constraints and can therefore be ruled out 
at the 99% C.L. 

Finally, we note here that M*/Lb and M tot / Lb{< 
-REinst) are correlated through Lb, which is only a scal- 
ing factor, useful for the physical interpretation, but ir- 
relevant for the lensing+dynamical analysis. Thus it is 
preferrable to express our results in terms of the fraction 
of dark matter /dm = (1 — M*/M to t), or equivalently in 
terms of the fraction of luminous matter /* = 1 — /dm- 
In the isotropic case, the range of dark-matter mass frac- 
tions inside the Einstein radius is /dm(< -REinst)=0.37- 
0.72. Translating this into a mass fraction inside R c is 
slightly model-dependent. However for 7=0-1, the range 
of dark-matter mass fractions is /dm(< -R e )=0.15-0.65, 
with at most a 10% change (both ways) in the value of 
/dm between 7=1 and 7 = 0. This confirms our conclu- 
sion that all five E/S0 lens galaxies at z w 0.5 — 1 have 
massive dark-matter halos, even well inside the luminous 
component. Note that these values for the dark-matter 
fraction arc significantly higher than the limits obtained 
from lensing statistics for adiabatically contracted lenses, 
/dm(< Re) < 33% (95% CL; Keeton 2001). 
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Fig. 9. — The cosmic evolution of the stellar mass-to-light ratio 
(for an isotropic velocity ellipsoid, rt = 00). The solid points are 
the results from the joint lensing & dynamics analysis, whereas the 
shaded region is the independent measurement via the Fundamental 
Plane from Treu et al. (2002). We emphasize that the latter is not 
based on lens galaxies and the former not on the FP measurements 
of M*/Lb from the five lenses as shown in Fig. 4. 



6.4. The evolution of M*/Lb from lensing & dynamics 

The sample of E/S0 lens galaxies reaches a large enough 
rcdshift to afford a direct measurement of the evolution of 
the stellar populations of E/S0 without including any con- 
straint from the FP. In Figure 9, we have plotted M*/Lb 
(Table. 6) as a function of redshift. The stellar mass-to- 
light ratios of E/S0 galaxies at z w 0.5-1.0 are significantly 
smaller than in the local Universe (on average 2.3 ± 0.6 
M q /Lb jQ versus 7.3 ± 2.1 M & /Lb,q) implying consider- 
able ageing of the stellar populations in the last 4-8 Gyrs. 



In Fig. 9, we also compare our direct lensing+dynamics 
measurement of the stellar mass-to-light ratio (Table 6) 
with the indirect measurement obtained from the evolution 
of the FP of field E/S0 (Treu et al. 2002). The agreement 
is very good, consistent with a scenario of pure luminosity 
evolution of E/S0 from z <~ 1 to today. This measurement 
rules out scenarios predicting strong evolution of the in- 
ternal structure of E/S0 galaxies with redshift, where the 
virial coefficient relating a 2 R e to the stellar mass would 
change significantly with cosmic time. 

If we use the value of (M*/Lb)'s from Table 6 to 
determine the evolution of the stellar mass-to-light ra- 
tio with redshift, we find that the average evolution is 
(dlog(M / Lb) / dz) = — 0.75±0.17, in good agreement with 
the results from Sect. 3.1, which gave d\og(M / Lq) / dz = 
— 0.72±0.10 based on the FP. A disagreement between the 
two independent results would have implied that cither the 
FP is not a good method to derive M*/Lb evolution or 
that our lensing+dynamical analysis is faulty. 
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Table 6 







Summary of lens/dynamical model results. 


II. 




Galaxy 


TV if 

M E 
(10 10 M Q ) 


(M tot /L B )<i ? , E 
M /Lb,0 


(M*/L B )™ (M,/L B )i s d 
M Q /L Bi0 M /L BjQ 


/dm ( 1S0 ) 
(< i?Einst) 


/dm ( 1S0 ) 

(< -REinst) 


0047 

C0302 

H1417 

H1543 

2016 


40.6±2.0 
67.0±6.7 
70.6±7.0 
3.4±0.7 
110.0±11.0 


5.4±0.5 
4.8±0.5 
5.0±0.5 
2.7±0.5 
8.0±0.8 
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Notes: The mass-to-light ratio is marginalized over 7 >0. 



but also in many application of lensing. For example, cos- 
mological parameters from lens statistics are generally ob- 
tained assuming that lenses are isothermal (e.g. Turner et 
al. 1984; Fukugita et al. 1990; Kochanek 1996; Helbig et 
al. 1999; Falco et al. 1999; Chae et al. 2002; Mitchell et al. 
2004). Similarly, the Hubble Constant from gravitational 
time-delays (Refsdal 1964) is typically obtained assuming 
isothermal mass density profile as well (e.g. Kundic et al. 
1997a; Schechter et al. 1997; Impey et al. 1998; Biggs et 
al. 1999; Koopmans & Fassnacht 1999; Koopmans 2001; 
Kochanek 2002; Wucknitz et al. 2003). Also estimating 
the central stellar velocity dispersion a from lens models 
depends on the assumed mass model. Hence, how justified 
is the isothermal approximation? 

1.5 



t — r 1 — 1 — 1 — 1 1 1 — r 



T~ 



aniso 



6.5. The inner slope of the dark-matter halos 

Since the values of M^/Lb agree between determina- 
tions from the FP and lensing plus dynamics, we can feel 
confident that their measurements can be combined, as 
shown by the solid lines in panels (a) and (b) in Fig. 7. 
The posterior probability distributions function of 7 - 
marginalized over M*/Ls, including the FP constraint - 
for the five lenses are shown as solid colored lines in Fig- 
ure 10. 

For individual lenses these posterior probability distri- 
bution functions imply upper limits on 7 between ~1 
and ^1.5, i.e. consistent with the inner cusps predicted 
by cosmological simulations, if the collapse of baryons to 
form stars did not significantly steepen the dark-matter 
halo (see TK02, KT02; see also Loeb & Peebles 2003; 
Sand et al. 2002, 2004; El-Zant et al. 2003; Nipoti et 
al. 2004). Somewhat tighter confidence limits are ob- 
tained when combining the measurements (dashed line): 
0.97 < 7 < 1.46 or most likely 7 = 1.3lo' 4 (isotropic) 
< 7 < 0.62 (anisotropic) at 68%CL and are 0.39 < 7 < 
1.59 (isotropic) < 7 < 1.26 (anisotropic) at 95% CL. In 
conclusion, the slope of the dark-matter halos is definitely 
flatter than isothermal, and ranges between the value pre- 
dicted by numerical simulations and zero, depending on 
anisotropy. Requiring consistency with numerical simula- 
tions implies that (i) significant radially anisotropic models 
(i.e. Ti w R e ) are ruled out and (ii) dark-matter halos do 
not significantly steepen during baryonic collapse. More 
stringent statements cannot be made at this stage, because 
uncertainties related to the orbital properties of the stars 
dominate the error budget 12 . 

7. THE HOMOGENEITY OF EARLY- TYPE GALAXIES 

Now that we have homogeneously analyzed a sample of 
five lenses at z ~ 0.5 — 1.0, we can start to look into the 
general properties of the E/S0 lens galaxies 13 . 

7.1. Are E/S0 galaxies isothermal? 

The first important question that we want to discuss is 
the average total mass density profile of lens galaxies, i.e. 
what is the distribution function of 7' (see Section 6.2). 
This is relevant not only in terms of formation scenarios, 

12 An additional source of uncertainty is due to the mass profile of the luminous component: the contours of 7' shift towards slightly lower 
values by adopting a Jaffe model for the luminous component. 

13 We note that the lens galaxies are not a statistical sample by any means, so the results should be interpreted as an exploration of the variety 
of possible behaviors rather than in a rigorous statistical sense. 




Fig. 10. — Posterior distribution functions for 7 for the individ- 
ual lenses and the joint probability. The case for isotropic orbits is 
shown in the lower panel, the case for radially anisotropic orbits is 
shown in the upper panel. 

The main conclusion from Section 6.2 is that there ap- 
pears to be intrinsic scatter in the values of 7' for lower 
values of i?Einst/^e- A possible trend is seen, with 7' in- 
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creasing with i?Einst/-Re and reaching 7' 2 at large rel- 
ative radii. However, based only on five systems, it is 
dangerous to interpret this as a real physical trend. For 
example, the two other lenses for which we have done a 
similar analysis PG1 115+080 (TK02b) and B1608+656 
(K03)- not selected as part of the LSD Survey 14 - have 
i > 2 (7' = 2.35 ± 0.1 ± 0.1 and 7' = 2.03 ± O.'l ± 0.1 re- 
spectively) . By including those in the plot (open triangles 
in Figure 8), the scatter in 7' would increase at low values 
of i?Einst/-Re, erasing any apparent trend. The straight 
average for this extended sample becomes (7') = 1.9 ±0.1 
with an rms of 0.3 for the isotropic case (r* = 00). Hence, 
with a larger sample in hand, one might indeed find that 
7' approaches w 2, but also that there is significant intrin- 
sic scatter (see also Rusin et al. 2003, who statistically 
build a mass profile from a sample of lens galaxies and 
find (7') = 2.07 ± 0.13). Is such a large scatter expected? 
At the radii of interests, typically of the order of the effec- 
tive radius, we can expect scatter for at least two reasons. 

First, E/S0 galaxies are typically found in or near groups 
or clusters and we expect that E/S0 lens galaxies will 
be in a similar environment (see discussion on the C0302 
lens model). A group dark-matter halo might be present 
around the lens (a cluster will be external, since they are 
critical themselves and would produce much larger image 
separations). If the inner mass slope of the group halo is 
shallower than isothermal (as typically found in the inner 
regions of groups/clusters; e.g. Ettori et al. 2002; Kelson 
et al. 2002; Gavazzi et al. 2003; Kneib et al. 2003; Sand et 
al. 2002, 2004), it could introduce a 'floor' of dark matter 
that will result in 7' < 2 (corresponding to the rising ve- 
locity dispersion profiles observed in some elliptical galax- 
ies at the center of clusters, e.g. Dressier 1979; Carter et 
al. 1985). Arguments suggesting that H1543, C0302 and 
possibly H1417 might be located in or near groups were 
presented in the previous sections. We emphasize that 
this is different from an intervening mass-sheet, since the 
group is physically centered on or close to the lens galaxy 
and thus also affects the stellar dynamics. On the other 
hand, if the Einstein radius is small enough, the total mass 
distribution will become more dominated by the luminous 
component, which is typically steeper than isothermal and 
results in 7' > 2. At face value, the trend seen in Figure 8 
is opposite to what we would expect: for small i?Einst/-R e , 
where baryon dominate we would expect 7' > 2, while at 
larger radii, where the group halo dominates, we would 
expect 7' < 2. Given the present size of the sample of 
five systems, the argument is inconclusive. The possible 
trend could simply result from C0302 and H1543 being 
in or near groups, and thus having a larger "dark-matter 
floor" then the other three lenses. However, this would 
be the opposite trend to e.g. PG1115+080 that shows a 
steeper than isothermal mass profile (TK02) and is also 
close to a group (Kundic et al. 1997). A possible mech- 
anism to explain steepening of the mass density profile of 
the lens through interaction with a group/cluster could be 
tidal stripping (e.g. Natarajan et al. 2002; TK03) of the 
outer halo. Thus, if the lens was located at the center of 
the group/cluster, where a significant amount of cluster 
dark matter is present the observed profile could be flat- 
ter than isothermal because of the "dark-matter floor". 



Viceversa, if the lens had been through the cluster/group 
center, deep enough to experience tidal stripping, but was 
observed on an outbound orbit far enough from the peri- 
center, the "dark-matter floor" would no longer be relevant 
and a steeper profile would be observed. 

Second, could the internal structure of luminous and 
dark matter in these E/S0 galaxies alone explain a trend 
or large scatter? In the local Universe, we know that 'ro- 
tation' curves of early-type galaxies show quite a variety 
of slopes within 1-2 effective radii (e.g. Bertin et al. 1994; 
Gerhard et al. 2001; Romanowsky et al. 2003), from rising 
(7' < 2) to declining (7' > 2). We might expect that the 
average and scatter of 7' would depend on the ratio be- 
tween the Einstein radius and the effective radius and the 
fraction of dark matter /dm contributing to the mass in- 
side the Einstein radius. The larger this ratio, the smaller 
the effect of baryons, and the more 7' will be a probe of 
the dark-matter effective slope (possibly including some 
effects of nearby clusters or groups). In this case 7' will 
increase in three stages (see Fig. 11): First, for r < r* (see 
Section 5) the slope is typically dominated by luminous 
mass and will have 7' < 2 for an i? 1 / 4 profile. Second, 
for r > r* a transition takes place where dark matter be- 
comes more prominent and the combination of luminous 
plus dark matter add to an effective slope of 7' w 2. How- 
ever, in the case that the fraction of dark matter (/dm) 
inside i?Einst is relatively small, one can expect that a very 
rapid transition can occur to 7' > 2 (e.g. PG1 115+080), 
if the luminous component remains dominant at r > r* 
for i?Einst somewhat larger than r*. Clearly this is very 
sensitive to /dm and i?Einst/^e and small variations in 
their value can induce large fluctuations in 7' for lenses 
with REinst/Re < 1. For larger values of i?Einst/-R e , the 
larger radial range covered will result in less scatter for 
the same changes in /dm- Third, around the break radius, 
a transition is expected from the region where 7' w 2 to 
a dark-matter dominated regime where 7' w 3 (see e.g. 
Seljak 2002 and Kneib et al. 2003 for a discussion of mass 
distribution at large radii from weak-lensing studies). 

It seems to us that most likely a combination of the ef- 
fects discussed above is required to interpret the observed 
trends and scatter in 7'. Group and cluster halos exist and 
must necessarily play a role (e.g. Lehar et al. 1997; Kundic 
et al. 1997a&b; Tonry 1998; Tonry & Kochanek 1999, 
2000; Blandford, Kundic & Surpi 2001; Keeton, Christlein 
& Zabludoff 2000; Fassnacht & Lubin 2002; Johnston et 
al. 2003). However, a model where galaxies are sim- 
ply isothermal and appear more shallow if embedded in 
a group/cluster halo, is not sufficient to explain the obser- 
vations. Some degree of internal scatter in the properties 
of the dark-matter halos of early-type galaxies is required. 
This scatter could be the result of complex and hierarchi- 
cal formation history and baryonic cooling history, and/or 
could be related to environmental effects (TK02b; Treu et 
al. 2003; Natarajan, Kneib, & Smail 2002). 

Ultimately, whatever the underlying cause or interpre- 
tation, we can not escape the conclusion that the inner 
total mass profile of E/S0 galaxies at z w 0.5 — 1.0 is on 
average close to or slightly more shallow than isothermal 
but also that there is a significant intrinsic r.m.s. scatter 
in 7' of up to ~0.3 (i.e. ~15% in density slope or ~30% in 



Because they fail to meet the criteria of favorable contrast between the lens and the source. 
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surface density slope). 

7.2. Lensing implications of inhomogeneity 

In the previous subsection we concluded that E/SO lens 
galaxies are close to isothermal on average, but not quite, 
and that there is a significant intrinsic scatter in the power- 
law slope of their total inner mass profiles. Here, we will 
briefly discuss the consequences of our finding on three 
important applications of gravitational lensing: (i) The 
determination of the Hubble constant from gravitational 
time delays; (ii) The determination of the star formation 
history of E/SO galaxies from image-separation estimates 
of the FP; (iii) The determination of the cosmological pa- 
rameters from lens statistics. 

Let us first consider the determination of the Hubble 
Constant from gravitational time delays and consider a 
lens with unknown mass profile, modeled as a singular 
isothermal ellipsoids. The observed rms scatter in 7' of 
~0.3 translates into a systematic uncertainty on Ho of 
~30%, in addition to other uncertainties. For example, 
an rms of ^20kms _1 Mpc -1 can be expected if the true 
value is, say, Ho=65kms _1 Mpc -1 . This range of 45-85 
kms -1 Mpc -1 covers the vast majority of determinations 
of Ho that assume isothermal mass profiles (e.g. Schechter 
et al. 1997; Impey et al. 1998; Biggs et al. 1999; Koop- 
mans & Fassnacht 1999; Kochanek 2002a&b; Koopmans 
et al. 2003; Wucknitz, Biggs & Browne 2003) and could 
therefore in principle explain the mutual inconsistency be- 
tween the inferred values of H . Different samples based 
on only a few lens systems could therefore lead to com- 
pletely different conclusions (e.g. Koopmans & Fassnacht 
1999; Kochanek 2002a&b), if 7' is not determined for each 
system directly. 

In TK02b and K03, we applied the lensing and dy- 
namics analysis described above to two lens systems with 
measured time delays, finding power-law slopes of 7' = 
2.35 ± 0.1 ± 0.1 for PG1115+080 and 7' = 2.03 ± 0.1 ± 0.1 
for B1608+656, respectively, which lead to values of H 
of 59^7 2 ± 3 kms" 1 Mpc" 1 and 75±£ ± 4 kms" 1 Mpc" 1 . 
In the case of PG1115+080, a 35% increase was found for 
the value of Ho from that expected from an isothermal 
model with 7' = 2 (i.e. Ho=44kms _1 Mpc -1 ; Impey et al. 
1998). A difference in slope of 0.3 between these systems 
is fully consistent with the r.m.s. scatter in 7' that we find 
in our sample and thus 'in hindsight' may not totally be 
unexpected. 

Indeed, in this paper we have presented three lens galax- 
ies with deviations of 7' in the opposite direction. Those 
systems would lead to severe overestimates of H if they 
were assumed to be isothermal. As we have stressed be- 
fore, one can overcome these problems by directly mea- 
suring the mass profile with a combination of lensing and 
stellar dynamics, or other methods (e.g. Wucknitz 2003). 
The statistical approach (i.e. lens galaxies are "on aver- 
age" isothermal) is not (yet) satisfactory, since the average 
value of 7' and its scatter are poorly determined at present. 
For example, if the current samples of about 4-5 lenses 
(e.g. Koopmans & Fassnacht 1999; Kochanek 2002a&b) 
were drawn from a Gaussian distribution of slopes with 
(7') = 2 and a l-cr width of 0.3, one would most likely 
find many 'outliers' (based on errors that do not incorpo- 
rate the uncertainty in the mass profile) with very low or 



high inferred values of Ho. The distribution of 7' for a 
sample of lenses might also depend on unknown selection 
functions and it appears therefore preferable to obtain a 
direct measurement of the mass slope (7') for lenses with 
time-delays. 

Let us now turn our attention to the star formation his- 
tory of early-type galaxies. Can we use multiple image sep- 
aration to estimate the central stellar velocity dispersion 
of E/SO galaxies to construct a Fundamental Plane of lens 
galaxies and measure the evolution of their mass-to-light 
ratio (Kochanek et al. 2000; Rusin et al. 2003; van de Ven 
et al. 2003)? What is the accuracy of this approximation? 
Based on our sample, we find that (/sie) < 1- Hence, if 
we had used the isothermal model to determine the central 
velocity dispersion of our lenses from image separation, we 
would have overestimated their effective mass (a 2 R e ) and 
underestimated the evolution of M^/Lb- If (/sie) was 
rcdshift independent, it would cancel out by fitting simul- 
taneously the local intercept (as done by R03), impacting 
the measurement only through the increased uncertainty 
due to the intrisic scatter in (/sie)- However, this could 
be a redshift dependent effect, since, for example, the ratio 
of the Einstein radius to the effective radius could depend 
on the redshift of the lens, and therefore it could mimic 
evolutionary trends. Hence, if the results from our sample 
of five E/SO lens galaxies hold for the larger sample of lens 
system, it could explain why most direct measurements 
(e.g. Treu et al. 2002; Gebhardt et al. 2003; van der Wei 
et al. 2004; this paper) indicate a slightly faster evolution 
of the Fundamental Plane of E/SO galaxies with redshift 
than those based on lens-estimates of a (e.g. Rusin et al. 
2003; van de Ven et al. 2003). If this difference - at the mo- 
ment only marginally significant - was confirmed by more 
precise measurements, it would be interesting to reverse 
the argument. If the lens-based estimate of the evolution 
of the FP is slower than the direct measurement, then 
the power-law slope of lens galaxies is on average flatter 
than isothermal and (/sie) < 1- The difference in the 
evolutionary rates would provide another measurement of 
(/sie) and more importantly on (7'). 

Finally, we shortly discuss the effects on lensing statis- 
tics. If galaxies have (/sie) < 1 and 7' = 2, then 
statistical models that use isothermal lens mass models 
with osie = o- will tend to underestimate the lensing 
cross-section of a population of lenses, since the correct 
csie = c//sie > o~. However, since we have found that 
/sie < 1 for many lenses might actually be associated with 
7' < 2, the effect is not so clear, since galaxies with a more 
shallow mass profiles have lower lens cross-sections for a 
fixed enclosed mass. On the other hand, they are more 
massive when normalized to the same a (i.e. a shallower 
profile lowers a for a fixed mass, hence to increase it to the 
same value as for a steeper profile, its mass needs to be 
increased). We are generally in the latter situation, since 
statistical models are typically normalized to an observed 
distribution function of a. It becomes even more complex, 
however, since all these effects need to be integrated over a 
distribution function of 7' and effects of the magnification 
bias as function of mass profile need to be accounted for. 
The full treatment of this problem goes beyond the scope 
of this paper and is left for future research. 

To summarize, the main conclusion of this section is 
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Fig. 11. — Density profiles of the five high-z LSD lens systems, representing the most likely isotropic (r; = oo) model for a dark-matter 
halo (blue line) with 7 = 1 (pdm <x r~~<) and a Hernquist stellar mass density profile (red line). The total mass density profile is plotted as a 
solid black line. The density is in units of 10 10 Mq kpc — 3 and the dashed line indicates the Einstein radius. The two slopes for 7' = 1 and 2 
are indicated for reference. 



that a simple isothermal model is probably not appropri- 
ate for precision measurements, especially when only small 
samples are available, as it is typically the case. The ob- 
served scatter does not seem surprising, given, for exam- 
ple, the wide range of velocity dispersion profiles observed 
for early-type galaxies, from declining to increasing with 
radius (see e.g., Gerhard et al. 2001; Kelson et al. 2002). 
Unfortunately the samples are still too small and the selec- 
tion functions too poorly characterized to find out what is 
the distribution of total mass density profiles and whether 
the distribution for lens samples differs from the one from 
morphologically or X-ray selected samples. As far as the 
lenses are concerned, at this stage the average slope ap- 
pears to be slightly shallower than isothermal, in marginal 
contrast with other determinations, using different meth- 
ods and samples (e.g. Rusin et al. 2003). Although the 
difference is relevant for precision measurements (e.g. of 
Ho) based on statisical arguments, we believe that what 
is remarkable is the relative agreement on the peak of the 
distribution being around 7' = 2, given the small size of 
the samples and the unknown sample selection biases. A 
larger number of lenses with precisely measured mass den- 
sity profiles is needed to make further progress on these 
issues. 

8. SUMMARY 

We have presented new spectroscopic measurements for 
three gravitational lens systems C0302 (z = 0.938), H1543 
(z = 0.497), and H1417 (z = 0.81) as part of the Lenses 
Structure & Dynamics (LSD) Survey. Long integrations 
with ESI at the Keck-II Telescope have yielded extended 
stellar velocity dispersion profiles of all three lens E/SO 
galaxies, extended approximately to the effective radius. 
A blue spectrum taken with LRIS-B has revealed the red- 
shift of the lensed arc in H1543 (z = 2.092). Together with 
two previously published systems, MG2016 (z = 1.004; 
TK02; see also Koopmans et al. 2002) and 0047 (z = 0.485; 
KT03), this paper presents the analysis of the current 



high-redshift sample (z w 0.5 — 1.0), consisting of five 
pressure-supported E/SO galaxies. 

The spectroscopic data have been combined with surface 
photometry from Hubble Space Telescope (HST) archival 
images to study the evolution of the stellar populations 
via the evolution of the intercept of the Fundamental 
Plane (FP). For the sample of five LSD lenses we find 
d\og{M/L B )/dz = -0.72 ±0.10, i.e. 1.80±0.25 magni- 
tudes of dimming between z = 1 and z = 0. In a pure 
luminosity evolution scenario, this measurement can be 
interpreted as the results of a relatively young luminosity- 
weighted age of the stellar populations. A scenario were 
most of the stars were formed at high redshift (> 2) while 
a small fraction of stars (~ 10%) is formed in secondary 
bursts between z = 1 and z = 0, provides a simple expla- 
nation for this result, as well as several independent pieces 
of evidence (evolution of the Luminosity Function of E/SO 
galaxies; spectroscopic evidence of recent minor episodes 
of star formation in distant E/SO; properties of local E/SO; 
sec the discussion in Section 3) 

New gravitational lens models of the three systems 
C0302, H1417 and H1543 have been presented. H1417 
has the classical 'quad' morphology and can be success- 
fully modeled with a Singular Isothermal Ellipsoid (SIE) 
mass distribution, with an external shear and a gradient in 
the local convergence that roughly aligns with the galaxy 
major axis and which we interpret as an internal asym- 
metry in the galaxy. In contrast, the lens systems C0302 
and H1543 are characterized by a source lensed into ex- 
tended arc-like features. These systems are modeled with 
an algorithm that allows for a non-parametric image re- 
construction (see Appendix), incorporating some of the 
techniques by Wallington et al. (1996) and Warren & 
Dye (2003). Both lenses are successfully modeled with 
a SIE mass model with external shear. In the case of 
H1543, a nearby (4'.'7) massive galaxy at the same redshift 
as the main lens is included in the lens model as a Singular 
Isothermal Sphere (SIS). The Einstein radii (i?Einst), SIE 
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velocity dispersions (ctsie) ano - enclosed masses (MEmst) 
of the three lenses are used to perform a joint lensing and 
dynamical analysis, with the following results. 

1. Constant mass-to- light ratio models (i.e. mass fol- 
lows light) are rejected at better than 99% CL for all 
five E/SO lens galaxies. A dark-matter halo with a 
mass density profile flatter than the luminous com- 
ponent is needed in all cases. The fraction of dark 
matter inside the Einstein Radius (/dm) is 37-72% 
(isotropic stellar orbits) and 15-65% inside the ef- 
fective radius. 

2. Modeling the total mass density profile of the lenses 
as a single power law density distribution p to t oc 
r~ 7 , the effective slope 7' is found to be on average 
somewhat smaller than isothermal, i.e. (7') = 1.75 
with and rms scatter of 0.20 (for isotropic veloc- 
ity ellipsoid; (7') = 1.57 with rms 0.35 for radial 
anisotropy) for our sample of five lenses. If we in- 
clude the two others systems that have a similar 
analysis, these values increase by ~0.15 and the rms 
increases to ^0.30. Consistent with these findings, 
the ratio /sie = er/csiE between central velocity 
dispersion and velocity dispersion of the Singular 
Isothermal Ellipsoid (SIE) mass model that best fits 
the lensing constraints is (/sie) = 0.87 with an rms 
scatter of 0.08. 

3. The average mass-to-light ratio of the luminous 
component (M*/Lb) = (2.3 ± 0.6) ft.65 M q /Lb, q is 
smaller than the average value for early-type galax- 
ies in the local Universe (7.3 ± 2.1) /i 65 M q /L_b, , 
consistent with passive evolution of a relatively old 
stellar population. The stellar mass to light ra- 
tio obtained from the joint lensing and dynamics 
analysis is found to evolve as (d\og(M / L&) / dz) = 
—0.75 ± 0.17 in excellent agreement with the inde- 
pendent measurement obtained via the Fundamen- 
tal Plane. 

4. The most precise constraints to date are obtained 
on the inner slope of the dark-matter halo 7. We 
find the following 68% confidence limits: 7 < 0.58 
(anisotropic velocity ellipsoid with = R e ) and 
0.93 < 7 < 1.48 or 7 = 1.3±g;| (68% CL) for n = 
00. Thus, our data are consistent with CDM numer- 
ical simulations (with 7=1-1.5) only if the veloc- 
ity ellipsoid is not significantly radially anisotropic 
and baryonic collapse, during galaxy formation, did 
not significantly steepen the mass density profile 
as would be expected in simple adiabatic contrac- 
tion scenarios (Blumenthal et al. 1986; Mo, Mao & 
White 1998; Keeton 2001; Kochanek 2002). 

9. CONCLUSIONS 

In conclusion, the following picture seems to be emerg- 
ing from the analysis of the complete high-redshift LSD 
sample. High redshift early-type galaxies are approxi- 
mately isothermal ellipsoids, but not exactly. Our cur- 
rent sample seems to indicate that on average the effective 
slope of the mass density profile inside the Einstein radius 
might be slightly smaller than 2, i.e. total mass density 



profile flatter than isothermal, albeit 7' = 2 is generally 
within the range of the distribution (e.g. the first two ob- 
jects we analyzed were almost exactly isothermal, TK02, 
KT03). A possible cause for departure from homogeneity 
could be the environment of early- type galaxies. Contri- 
butions from a relatively flat group or cluster dark-matter 
halo could introduce a "floor" of mass at the position of the 
lens causing the total mass density profile to appear effec- 
tively flatter. Independent external evidence (such as the 
presence of nearby galaxies at the same redshift) indicates 
that possibly all three lenses for which we found 7' < 2 
might be members of a group, and therefore this mecha- 
nism would appear to be a likely explanation, at least for 
this sample. However, this mechanisms does not appear 
to be sufficient to account for all the observed scatter, be- 
cause most of the scatter is observed for i?Einst/-R e < 1, 
where the galaxy mass distribution should dominate over a 
group halo, and because there are cases where 7' > 2 is ob- 
served (such as PCI 115+080 and B1608+656; the former 
also has a nearby compact group). We therefore conclude 
that the scatter in the total mass density profile is asso- 
ciated with intrinsic scatter in the ratio of dark matter 
to luminous matter in the inner regions of high redshift 
early- type galaxies, similar to what is observed in the lo- 
cal Universe (Bertin & Stiavelli 1993; Bertin ct al. 1994; 
Gerhard et al. 2001). 

Regardless of the physical interpretation, the observed 
scatter in the effective slope (i.e. 10-15% in density slope) 
still implies a remarkable degree of structural homogene- 
ity between early-type galaxies from a galaxy formation 
point of view. Whatever the details of the mass assembly 
and star formation history, E/S0 galaxies end up being 
close to isothermal, for example expected in (incomplete) 
violent relaxation scenarios (Lynden-Bell 1967; Shu 1978; 
van Albada 1982; Bertin & Stiavelli 1993; Hjorth & Mad- 
sen 1991; 1995). From an evolutionary point of view, the 
intrinsic scatter in the mass density profiles of high red- 
shift E/S0 galaxies does not seem to be much different 
from that of local E/S0 galaxies, providing no evidence 
for much structural evolution within the last 4-8 Gyrs. 

Although the homogeneity is remarkable from a galaxy- 
formation point of view, the observed scatter in 7' is large 
enough that the isothermal approximation might not be 
good for some applications, particularly when they depend 
critically on the mass slope. Meaningful examples include 
the determination of Ho from gravitational time-delays, 
where a scatter of ^0.3 in 7' translates into a scatter of 
~30% in H for a given time delay (Saha 2000; Wucknitz 
2002; Kochanek 2002; TK02b). Thus, it appears necessary 
to use external information - such as internal kinematics 
- to pinpoint 7' and Hq to a level of accuracy (10-15%) 
competitive with other methods (Koopmans et al. 2003; 
see also Wucknitz et al. 2003). A precision measurement 
based on statistical assumptions of 7' will probably have 
to wait until the distribution of 7' is observationally well 
characterized, and a large enough sample of lenses with 
time delays is available to minimize variance. 

The other main result of this paper is the decomposition 
of the total mass distribution into a luminous component 
and a dark-matter halo. High-redshift early-type galax- 
ies are inconsistent with constant mass-to-light ratio mass 
models. The mass-to-light ratio has to increase signifi- 
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cantly with radius, consistent with the presence of dark- 
matter halos flatter than the luminous component (dark 
matter makes up a substantial fraction of the mass of the 
lenses inside the Einstein Radius, of order 40-70%). 

The mass-to-light ratio of the luminous component is 
smaller than in the local Universe, consistent with pas- 
sive evolution of a relatively old stellar population. The 
agreement between the evolution of the stellar mass to 
light ratio measured directly and that determined from 
the evolution of the intercept of the Fundamental Plane is 
another argument against strong structural and dynamical 
evolution. If early-type galaxies had changed their mass 
distribution significantly between z ~ 1 and today, they 
would have to be doing so while preserving the mapping 
between velocity dispersion, radius and stellar mass. A 
simpler - and to our eyes preferable - explanation would 
be one involving no or little dynamical evolution. 

The precise mass density profile of the dark-matter ha- 
los is harder to constrain, since most of the mass at small 
radii is luminous. Nevertheless, our sample of five lenses 
allows us for the first time to set limits on the inner dark- 
matter density slope 7 that are interesting to compare with 
cosmological simulations. Our measurement can be recon- 
ciled with numerical cosmological only if (i) the velocity 
ellipsoid is not significantly radially anisotropic; and (ii) 
the collapse of baryons to form the galaxy did not steepen 
the dark-matter halo more than a few tenths in 7. The lat- 
ter constraint is clearly inconsistent with simple adiabatic 
collapse models (Blumenthal ct al. 1986) and suggests that 
different mechanisms are involved in the accretion of stars 
of the centers of the halos of early-type galaxies. This is 
further supported by lcnsing statistics results, which im- 
ply a low value of / DM (< R e ) < 0.33 (95% CL; Keeton 
2001). Assuming adiabatic contraction implies that the 
initial slope (7$) of the dark-matter halo (i.e. before con- 
traction) was shallower than predicted from ACDM mod- 
els (i.e. ji < 1) and introduces an additional inconsistency. 
Mechanisms such as those proposed for clusters (El-Zant et 



al. 2003; Nipoti et al. 2003), where stars form first in satel- 
lites, which are then accreted as dissipationless particles 
could provide the desired effect (see also Loeb & Peebles 
2003), although a detailed comparison of data with theory 
will have to wait until cosmological simulation including a 
realistic treatment of star formation become available. 
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APPENDIX 
A. DELENSING OF IMAGES ON A GRID 

Suppose we have a lensed image, d, on a grid (e.g. a CCD image) 15 in which all emission not associated with the lensed 
source has been masked and/or subtracted. Hence the grid should be a noisy, blurred and lensed representation of the 
true underlying source brightness distribution, which it is our aim to reconstruct. Suppose also that we can construct 
a lens-operator L (depending on the parameters of our lens model) which acting on a source grid, s, produces a lensed 
image of the source. Suppose further that a blurring operator, B, exists which acting on Ls produces a blurred lensed 
image. Putting this together, we have 

BLs = d t +n, (Al) 

where n represents the noise in the observed image and dt = d — n is the noise-free lensed image. Note that the size 
or shape of the image and source grids and their pixel sizes are irrelevant to the this problem. Furthermore, neither 
the grid nor the pixels have have to be rectangular and connected (e.g. the image grid could have a gap). As long as 
the value of each observed pixel can be written as a linear combination of source pixel values, the above equation is 
applicable. For simplicity, however, we will assume rectangular grids and pixels for the remainder of this paper. Naively, 
one would think that the solution to this problem, i.e. the source brightness distribution, can easily be found by inverting 
the above equation through § = (BL) -1 d. However, this is a notoriously ill-posed problem and noise in the observed 
image will typically lead to unacceptably poor reconstructions. We will not further discuss this solution. The problem can 
be regularized, however, suppressing the effects of n on the final source reconstruction (Thikonov 1965). Mathematically, 
we would like to find the source grid s and the parameters of the lens model (i.e. L) that minimize the following equation 

C(A) = ||BLS-d||| + A||HS|||, (A2) 
where H is a rcgularization operator and A determines the weight given to the regularization term. For simplicity we 
have also written (BL)/n — > (BL) and d/n — > d (see Press et al. 1992). The first term in Eq.(2) is simply the \ 2 term, 
whereas the second term regulates the "smoothness" of the final solution. In particular, if H = I, then the regularization 
term is simply the sum over the squared pixel values in the source grid. The latter term, however, can also regulate the 
smoothness of the derivatives in the source grid, or be replaced by a maximum likelihood or maximum entropy term (e.g. 
Whyth et al. 2003). The solution to Eq.(2) is the solution to the set of linear equations 

[(BL) T (BL) + AH] §= (BL) T d. (A3) 

This equation has a unique solution, thanks to the regularization term, and can be solved through standard techniques 
and using freely available linear algebra packages for large sparse matrices. 

B. THE LENSING & BLURRING OPERATORS 

The main problem that is faced in constructing a lens operator is the fact that an image pixel, when projected on the 
source grid, in general will not exactly coincide with a source pixel. One way to solve this (Warren & Dye 2003) is to 
determine the lensed image for a given source pixel of unit flux, blurred by the PSF and regridded to the image grid. 
This is repeated for all source pixels. The source is then reconstructed by finding the set of source-pixel weights (i.e. 
fluxes for each source pixel) and the corresponding lensed images, for which their linear combination best reconstructs 
the observed image. This requires inverting the lens equation and is numerically expensive for complex models. It also 
requires finding all lensed images for each pixel. There are no simple algorithms that guarantee this for complex mass 
models. However, one can also construct a lens operator that does not require the lens inversion. Before describing this, we 
first introduce some definitions. We assume that the source grid has K x L pixels (K columns and L rows) and similarly, 
the (as of yet unblurred) image grid, d', has M x N pixels. The values of each pixel are Sk,i and d' mn , respectively. Hence 
s = {s i= k+(i-i)K;} with i = 1...KL with k = 1...K and I = 1...L. Similarly, d' = Ls = {d'j =m+ i n _ 1 \ M } with j = 1...MN 
with m = 1...M and n = 1...N. It is now very easy to contract a particular implementation of the lens operator, L, 
which is a matrix of size KL x MN and entries kj. We emphasize that our choice is not unique, but it is simple and 
fast. First, each pixel j = m + (n — \)M is cast back on the source plane to a position yj, using the lens equation. If 
the position is outside the (pre-defined) grid, one simply continues to the next pixel. In general, however, the lensed 
image grid and source grid will be defined to overlap as much as possible to reduce redundancy, although this is not a 
requirement. Second, one determines the four pixels in the source grid that enclose yj, say (r + /i, s + v) for fi = 0, 1 and 
v = 0,1. Suppose that the position of (r, s) is y riS and (t, u) — yj — y rtS . The flux of pixel j can then be written as a 
linear combination of the four source pixel fluxes 

l l 

d '] = X! X] W ^< u S i=(r+t*) + (s+v-l)K (Bl) 
jt=0 u=0 

15 All two-dimensional grids are represented as vectors in which consequativc rows are placed behind each other. 



24 



Treu & Koopmans 



where 

w 0: o = (1 — i)(l — u) 
Wlfi = t(l - u) 
wo.i = tu 

w 1A = (1 - t)u. (B2) 

This is simply a bilinear interpolation of the four source pixel fluxes, but more complicated linear schemes can be 
constructed, although they will generate more entries in the lens operator (except if ones chooses to use an enclosing 
triangle). The four entries in the lens operator, at kj, for image pixel j are then the values of at i = (r+/j,)+(s+v—l)K 
for /i = 0, 1 and v — 0, 1. Finally, L contains at most AMN entries. 16 

The next step is to construct the blurring operator, B, that acts on the lensed image d'. As illustration, suppose that 
the PSF is a square grid of size (2H + 1) x (2H + 1) pixels with values with /i = —H...H and v = —H...H and peaks 
at p ,o- The sum of the PSF pixels adds to unity. The entries in the blurring operator, b g ^ are then simply the values of 
V^.y at g = {h + //) + (h + v - 1)M, if and only if 1 < (h + fi) < M and 1 < (h + v) < N for each h = 1...MN. Notice here 
that this method allows one to define a color-dependent blurring operator. In that case, for each h (i.e. pixel in the image 
plane), one can use a PSF with values p^^.h that depend on the local color of the pixel. We also note that an extinction 
correction can be done by setting the integral over the PSF, E M £„p^„ = e~ T < 1, where t is the optical depth due to 
extinction. (Note that this requires a color and extinction model, since both effects occur before blurring.) 

C. PRACTICALITIES IN THE OPTIMIZATION 

The pixel size in the source plane is set roughly by the largest pixel magnification in the image plane. If the source- 
plane pixels are chosen too large, the resulting image shows the effects of the mapping of individual source-plane pixels. 
If they are chosen too small, however, the source breaks up in "strings" of pixels that map closely to the image-plane 
pixels, but no flux in between (e.g. for H=I) Mathematically, these solutions are equivalent, but physically clearly not. 
We therefore set the pixel size such that the source does not tend to break up and adjacent image-plane pixels roughly 
map onto adjacent source-plane pixels. We then set A to an initially large value of typically ~0.1 and minimize C(A) by 
varying the lens- model parameters (see WD03 for details). We then lower the value of A slowly, continuing to optimize 
the model parameters, until a reduced x 2 of unity is reached. In general we find that the resulting mass-model parameters 
are relatively robust against changes in cither the source-plane pixel size and the value of A, as long as the source remains 
compact, and we refer to WD03 for a more thorough discussion of different choices. 



16 Note that L contains a fraction of < 4/(KL) non-zero elements. For a 100x100 image grid, this is <0.04%. It is obvious that sparse-matrix 
packages are required, or otherwise L alone would require ~1 GB of computer memory to store. 



